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ABSTRACT 


Tlie nonlinear equations of motion for a rigid rotor rastraine 
by three flexible springs representing, respectively, the 'flapping 
lagging and feathering motions are derived uci.u>v Lagrange’s equa- 
tions, for arbitrary angular rotations. ThcHa a: e redn.cc:-: to a 
consistent set of nonlinear equations using ncal t-.erris up to 

third order. The complete analysis is divided into tl-.ree parts 
A, B and C. 

Part A consists of forced response of tvv.i-dogree iiapping- 
lagging rotor under the excitation of pure greviraj-ir^n,:! field 
(i.e., no aerodynamic forces). Both forced osc::l‘! j. at i oiiS an we].! 
as parametric resonance are investig^lted u:->incj tr-e T-salanc 

method and solving the resulting nonlinear alcebraic oqua r iojis 
numerically by Nev/ton-Raphson iterative r.eohn i.que . i f f ecti of 
initial coning angle and flapping to lagging ii'.-og'.iHnoy ratio £ire 
discussed. For relatively small initial coning aiigle (about 9*^) 
the nonlinearity becomes softening spring type- and large coupled 
responses are possible for rotational frequencies signif ie-intly 
lower than the lagging frequency. 

In Part B, the effect of aerodynamic forces on tlie dyncunic 
response of two-degree flapping-lagging rotor is investigated.. 
Significant aerodynamic effects are found for some of the previous 
forced oscillations and parametric resonemces. Also, self-excited 
aerodynamic flutter instabilities are obtained after neglecting 
tlie gravity forces. Effects of vaious parameters like Lock number 
inflov.' ratio, initial coning angle, structural damping, etc. are 
discussed. .Also, the effect of a v;ind shear velocity gradient is 
investigated, and is found to produce little effect on the lagging 
response but appreciable effect on the flapping amplitude. 



In Part C, the effect of third degree of motion, feathering 
is considered. First, the forced response of flapping-lagging- 
feathering rotor under gravitational field and with wind shear 
flow is studied. It is found that even for relatively torsionally 
stiff rotor, the flapping amplitude is increased and the feathering 
response is appreciable. For the self-excited aerodynamic flutter 
instability, it is found that the feathering motion can reduce the 
linear instability speed appreciably. Also, the limit cycle flutter 
solution of a typical configuration shows a substantial nonlinear 
softening spring behavior. This reveals the possibility of sus- 
tained limit cycle flutter oscillations occurring well belov; the 
linefir instability speed if large enough disturbances arc. given to 
the rotor. 
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SECTION 1 
INTRODUCTION 


For centuries, wind energy systems have been used as sources 
of power in different forms like sailing of ships, pumping of 
water, grinding of grain, generation of electricity, etc. In 
the middle of the twentieth century, the interest in these sys- 
tems declined because of the poor economic viability with other 
power generation systems. However, with the recent energy crisis 
and also because of readily available technology of fixed and 
rotary v/ing aircrafts^ wind power is being considered as one of 
the potential sources of clean nondepleting energy. 

Many types of wind power systems have been tried in the past 
to trap the thinly distributed wind energy. Most of the wind 
machines can be broadly classified into three categories depend- 
ing upon the orientation of the axis of rotation, as Horizontal- 
axis wind rotor. Vertical-axis wind rotor and Cross-wind rotor. 

The Majority of these machines are based on the principle of 
Airfoil lift. At the present time, more attention has been given 
to the Horizontal-axis wind turbine partially because of its 
better theoretical understanding. So all the latter discussion 
and analysis is related mainly to tliis type of machine. (For 
general material see Refs. 1-8.) 

Because of the low density of air, large amounts of air must 
be tapped to provide an appreciable amount of power. It is 
generally true that the cost of pov/er produced is reduced with 
increasing size of wind driven plant, i.e., increasing the size 
of the rotor. However, one of the prime design problems of the 
big wind turbine system is the dynamics of rotor blades and the 
supporting tower and this has direct bearing on the operating life 
of the rotor. To keep low the cost of power production it is 
essential that the rotor should have long efficient life, i.e., 
subjected to less vibratory fatigue loads. Thus, there is need to 
understand the dynamic characteristics of the wind turbine system. 
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In some aspects, the dynamics of the big v;ind turbine is 
quite similar to that of the rotary wing aircraft. In the past 
three decades, a lot of research has undergone to understand the 
various dynamic or aeroelastic problems of helicopter and tilt 
rotor aircrafts, see for example Loewy's review paper [9] and 
other Refs. 10-23. A good deal of techniques developed in the 
formulation and analysis of the aeroelastic problems of rotary 
wing aircrafts can be used to study the dynamics of wind turbine 
rotors. However, all the results of rotary wing aircraft cannot 
be transformed directly for v/ind turbine because of differences 
in some of the parameters like rotational speed, tip velocity 
ratio, stiffnesses and v/eight properties, etc. In addition there 
are certain specific aspects of wind turbine dynamics which have 
to be looked into individually. For example, forced response of 
wind turbine blade under periodic forces due to gravitational 
field and sheared flov; effect and also impulsive forces due to 
tov/er shadov/ effect are quite important problems, particularly 
for big wind turbines. A good general picture of various dynamic 
problems concerning wind turbines is given in Refs. 24-26. 

There is little literature available related directly to 
the dynamics of v/ind turbine. Orraiston [27] has made a simple 
linear analysis from the uncoupled flap and lag equations for the 
forced response of v/ind turbine rotor under the excitation of 
gravity forcing function and also due to velocity gradient effect. 
The influence of blade number and hub articulation on the blade 
and tower stresses is examined and also the basic scaling relation- 
ships with respect to the length of the blade are discussed. Kaza 
and Hammond [28] has formulated the general linear flap-lag equa- 
tions for flutter stability applicable both to the wind turbine 
rotor v/ith velocity gradient as v/ell as helicopter rotor in for- 
ward flight. Two types of hinge sequences for flap-lag motions 
are used and the equations v/ith the periodic functions are solved 
using the Floquet-Liapunov method as well as the approximate 
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method (time averaging of periodic functions) . It was seen that 
the velocity gradient has little effect on the flutter boundary 
v;here as hinge sequence for flap-lag motions has a strong influence 
on the flutter stability of this two degree of freedom system. 
Friedmann [29] has derived the general coupled nonlinear flap- 
lag-torsional equations of motion for moderately large deflec- 
tions of a pretwisted cantilevered wind turbine blade with the 
incoming wind having velocity gradient as well as gust components 
in all the three directions. The methods to solve these equations 
are mentioned. Miller [30] has obtained the linearized version 
of the nonlinear flapping-lagging-feathering flutter equations of 
rotor by considering the motion to be small perturbations about 
possibly large static solution. The importance of various physi- 
cal quantities involved in the flutter and divergence of windmill 
blade is discussed. The effective damping plots are obtained for 
various configurations from the eigen analysis of the flutter 
equations. Dugundji [31] has given a good review of the v/hirl 
stability problem of wind turbine rotor mounted on a flexible 
tower. The general linear coupled equations of motions are 
derived for flapping-lagging rotor with two degrees of motion of 
tower head. The solution of these equations containing periodic 
coefficients using Floquet theory, for two bladed rotor particular- 
ly, is discussed. Some experimental results of small windmill 
model are given. In Refs. (26,32,33), the authors discuss the 
various aspects like design, fabrication, analysis, testing, etc. 
of 100 KW NASA Wind Turbine also discuss some of the dynamic 
problems pertianing to this wind turbine. 

For most of the aeroelastic analysis of rotors, the basic- 
ally nonlinear equations of motion are linearized by retaining 
only important static terms. Then it becomes much easier to work 
on the linearly coupled equations. However, there are some non- 
linear analyses in the literature. Young [34] has made a quali- 
tative analysis of the second order nonlinear equations of flapping- 
lagging rotor by the approximate method. Hohenerascr and Heaton [35] 
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have used the stepv/ise mamerical integration scheme to solve the 
second order flap-lag equations. Tong and Friedmann [14] has made 
an exhaustive nonlinear analysis of flap-lag as well as Flap-Lag- 
Feathering rotor by the multiple time scales perturbation method. 
Another method. Harmonic Balancing, is quite widely used in the 
linear dynamic analysis of the rotor mainly because of its simpli- 
city (e.g.. Refs. 15,22,36). Dugundji, etc. [37,38] have used 
Harmonic Balance Method to solve the nonlinear panel flutter equa- 
tions as well as to obtain the nonlinear forced oscillations 
response of the beams. 

In the present report, nonlinear dynamic analysis is made for 
an isolated blade of v;ind turbine with no tower interaction. The 
blade is assumed to be completely rigid and is restrained by three 
flexible springs at the hinge point representing, respectively, 
the flapping, lagging and feathering degrees of motion. It is 
further assumed that the blade c.g. , aerodynamic center and 
elastic axis lies at the quarter chord point and there is no 
variation of any of these along the blade axis. A particular 
hinge sequence of feathering first (from rotation axis) , flapping 
second, and lagging motion last is followed. However, one can 
expect different results with changed hinge sequence [17,28]. The 
equations of motions are derived using the energy approach (i.e. 
Lagrange's equations). Keeping nonlinearity up to third order, 
the consistent nonlinear differential equations are obtained. 

The complete analysis is divided into three parts. A, B and C. 

Part A consists of forced response of flapping-lagging rotor 
under the excitation of pure gravitational field. No aerodynamic 
forces are considered here. The blade can, however, have initial 
feathering angle setting. Both forced oscillations as well as 
parametric resonance are investigated. The forced oscillations 
response takes place at the frequency of the forcing function 
(i.e. rotational freq.) where as for parametric resonance the 
response frequency is one half the forcing frequency.. First, 
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simple linear solutions are worked out from the uncoupled flap 
and lag equation to get some basic understanding of the possible 
response of the blade. Then nonlinear limit cycle solutions are 
obtained for the flapping-lagging equations by applying the 
Harmonic Balance method and solving the resulting nonlinear 
algebraic equations nxamerically by Newton-Raphson iterative 
technique. These solutions are checked for stability to see 
whether they are physically existent or not. The stability check 
is snade by giving small perturbations to these steady solutions 
and studying the growth rate of these disturbances with time under 
the assumption of slowly changing functions. If the perturbations 
grow v;ith time means solution is unstable. The effect of initial 
coning angle and flapping to lagging frequency ratio on both 
forced response as well parametric resonance is investigated. 

The comparison of linear and nonlinear solutions near and away 
from resonance conditions is discussed. 

In Part B, the effect of aerodynamic forces on the two-degree 
flapping-lagging rotor is investigated. Quasi-steady airfoil 
theory is used to obtain the aerodynamic forces. First forced 
response of rotor is studied under the excitation of gravitational 
forcing field and in the presence of aerodynamic forces. Again, 
the nonlinear analysis of flapping and lagging equations is made 
like Part A for both forced oscillations as well as parametric 
resonance. The effect of various parameters like Lock number, 
inflow ratio, coning angle, structural damping, flapping to 
lagging frequencies ratio etc. on response amplitude is investi- 
gated. Then the self-excited flutter response of this torsionally 
rigid rotor is studied in the absence of gravitational forces. 

The equations of motion are the same as the first case except 
that all periodic terms are absent in these equations because of 
neglecting gravity forces. First simple linear analysis is made 
and then more rigorous nonlinear solutions are obtained by using 
the Harmonic Balance method. The nonlinear flutter solution is 
slightly different from that of forced response, here, for a 
knov;n lagging amplitude the solution is worked to obtain the 


5 



corresponding flapping amplitude, the flutter frequency and the 
stiffness of the configuration. The results are presented in 
the form of stability envelopes. The effect of inflow ratio, 

Lock member, coning angle, structural damping, hinge offset, etc. 
on the ciritcal flutter boundary are discussed. The behavior 
of limit cycle flutter amplitude with changing rotational speeds 
is also studied. In the end of Part B the effect of sheared 
flow on the forced response of flapping-lagging rotor is investi- 
gated. Expanding the velocity profile power law relation and 
retaining terms up to the second order, and comparing the elemental 
thrusts obtained from momentum theory and blade element theory, 
the inflow at any point'' is expressed in terms of inflow at the 
hub, blade azimuthal angle and the radial distance of the point. 

The equations of motion here get modified and these contain 
periodic aerodynamic terms. Again, by nonlinear analysis the 
effect of velocity gradient on the forced response of the blade 
is studied v/ith and v/ithout the gravity forces. 

In Part C, the effect of third degree of motion, feathering, ■ 
normal to the axis of rotation is considered, thus making the 
rotor a three degree of freedom system. The general equations of 
motion for this flapping-lagging-feathering rotor in the presence 
of gravity forces and with sheared flow are worked out. First 
simple solutions are obtained, then nonlinear limit cycle solutions 
are obtained for the forced response of the rotor by the Harmonic 
Balance Method. The nonlinear solutions are again checked for 
their well-posedness . The response amplitudes of a typical rotor 
configuration v/ith the three degrees of motions, for tvjo cases 
of with and without sheared flow, are compared with those of the 
same rotor with the feathering degree of motion locked. Then the 
self-excited flutter solutions for flapping-lagging-feathering 
rotor are investigated after neglecting gravity forces and also 
considering the uniform inflow. First the linear flutter analysis 
is made by assuming the motion to- be small perturbations about 
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some possibly large static solution. The damping plots are 
obtained for various configurations from the eigen analysis of 
the linearly coupled equations. Then nonlinear flutter analysis 
is made by the Harmonic Balance method. The behavior of the 
limit cycle flutter amplitude with changing rotational speeds is 
studied. 
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PART A; NO AERODYNAMIC FORCES 


SECTION 2 

NONLINEAR EQUATIONS OF MOTION 

The rotor blade v;ill be considered rigid with root hinges as 
shov/n below. The flapping and lagging hinges have the same offset, 
'e', and the C.G. of each blade cross-section is assumed to lie on 
the longitudinal ^ axis of the blade. No aerodynamic forces will 
be considred at this time, since the main purpose here v/ill be to 
assess the effects of the gravity forces on the rotating windmill 
rotor blade. 


Hinge Joint 


Rotation 

axis 



- 

- 0 

- 3 

- 


For convenience in setting up the nonlinear equations of 
motion, the no-offset-case (e=0) will be derived first, then the 
effect of the offset, e, will be added later. 

All equations and all subsequent calculations in this section 
v/ill apply for the hinge sequence shown above, i.e., feathering 0 
first, flapping B second, and lagging <}> last. 
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2 . 1 No Hinge Offset Present 


When no offset is present, the absolute location x, y, z of 
any point on the blade, n# can be defined in terms of four 
axis rotations involving the Euler angles ip, 8, B, and ip respectively. 
These are shown in the sketch below. 
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Multiplying out the various rotation matrices above, gives the 
following relation between the fixed axes x, y, 2 and the blade 
axes 5, n, ^ namely. 


I- 
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It should be noted that the inverse of the above square matrix 
is equal to its transpose. 

From the above relationship, the absolute velocities x, y, z 
of any point ri/ 5/ can be found by differentiation to be. 
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The kinetic energy T of this rotating blade is, 


T = i } dw (5) 

The evaluation of Eq. (5) is extremely tedious as it involves 
many trigonometric terms in the squaring process which later 
combine together. To circumvent this, an alternative expression 
is used for kinetic energy T based on the fact that the blade is 
in pure rotation about the origin. It can then be shown that the 
kinetic energy T is also. 


where u^, oj^, are the angular velocities of the blade about 

the Cf Hi C axes respectively. For a flat blade v/ith the section 

C.G. on the ^ axis, the product of inertia terms I,. Xl ^ 1 ^ r 0 , 

^ Cn nc 

hence only the first three terms need be retained. 

To evaluate the angular velocities u)_, u) , w , in terms of the 

.... n ? 

coordinate velocities i</, 0, B, <f>, one notes that the total vector 
velocity w can be written in tv/o ways. 
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From the rotation transformations between the various axes, 
the following relations exist between the unit vectors, 
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Placing these into 

(8), reducing all unit vectors to 4^, 

C t y t 

n c 


and comparing with 

(7) gives the angular velocities as. 
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where ij; = 57. is the constant rotation speed of the rotor. Using 
the above expressions for angular velocities, the kinetic enex-gy 
(6) for this blade with I_ “Z 1 ^ 'Z I - Z o becomes, 

nc CC 
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The internal potential energy U of this rotor blade arises from 
torsional springs kg, k^, k^ which are placed at the three hinges > 
of the hub. This gives, 

i- 4 .^ ( 4 ^- 4 s.) 
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where 0 . 3_/ <{>_ tire the initial settings for no spring moments. 

w ^ ^ 

The gravity forces give rise to an incremental work 6W as shown 
in the sketch below. 


A 



7 



ity ad’s- % 

I 

S'W = J 

m = V^aS5>/i>i 


axis 

(13) 


♦ 


The incremental displacement 6x can be expressed in terms of the 
incremental variables 60, 68, 6<}) using (2) as. 


— X St 


^ j5tfi 0 CtfvL 0 ^ S 0 (14) 

p C(A- ^ ^ ^ 

/iAyiO /Xvv ^ Cj^ G ^ ^ 

[ ^ , 


etc. 


T 
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since the blade section C.G. is assumed to lie on the C axis, 
then n = C = <? in (14) . Placing (14) into (13) and integrating 
gives 

I 

- G/q^o S'p -t ^15) 


where. 


Qa = 
- 

Q, - 


^ ^ ^ ^ CG^ 0 ^ CG^ (J? t 0 A^i/vv. 

^aX^B — 6ro.^ (^aJ^ ^ 

+ Ccr^0 

Co^ ^ ^ OerO, ^ ,4 aw <1?^ 


(16) 


X. 

=: t>tatic VvVi balance ==. I W\ ^ 
aboul" ^ o 


If aerodynamic forces were also present, additional aerodynamic 
moments Q^, Q^, would be added to (15). 

Gathering together the T, U, from (11) , (12) , (16) , and 

placing into Lagrange's equations 


f:^) - 2]1 + ih: = Q. 


(17) 
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N 


gives, after some algebra, the nonlinear large deflection equations 
of motion as. 


0 £g}M.a: 


e 


^ I ^ ^ ccrV- ^ t X •j p ^ 

** / \ ** 

^ \ ~ ) X ^ /<uax, z ~ 4 * X ^ ^ 

* * / 1 - 2 . \ 

“ 6p ^ CcrU. 4- 1^ /4 -<a<^ X-j y Z ^ 

0 4 (3 Z4 + p 4 ) ortiz4 ~X3.j|c^^ 

*” ^ (X-j," l^*) (^ aav.'Z^’ -t SX ^ j A^0/<iXvv ^ >i1a^ 

[T| cj? i- c<r=^0 + ^X-^ Ur^ c{> + X>^ /u^cf> -'1,-^c^B Co-^Z^J 

C-(T^0 A^uvuZ^ /4 XvuZ 4* ~ A^v 0 ^ c<jra. Z^ ■+ ^'• 

X ^ ^ "'* 4 ^ a^Z 9 p — (^X^/<u<w 4 ■*' X 7 j<i ^43 ^< 5 XtZ 0 

+ (^X^- I^) co<tZ0 Xwv^ -V I 3 ^Ze j 4- J?'!.0(e'-0s) 

~” ^ ^ C(5^ 0 ^ C^5<1.c|7 4- /li>y,<^^ 


(18) 
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^ Ec\UCL I 


- p4>(i,-ii) A^z<t> + et [(r^-rj')c^f «<2<)> + x^c^p] 

-SI6[( 

^ ^ B" Aa/vv- ^ xaA-vn. 7. 4" ^ 

■i- SZ^ + crsL Z4 c/j^e j 


1 

2. 


4- 51 


^ 1^— ^ 6-^ra- 4 X->j 0 A'^T. ^ "*■ ( A'^aZy 


4- X-« C^ 0 




Zp] Ap(p-p5) 

= ^-s j'-<Xv-'}' (>1^S £<«.<(’) — Cera.'f ( /ai-^p ««-t) 


(19) 


E 


UCL ^ 


^ ~ *0 - 0p ( I^-I^) c^p c^Z4 + T^ 

0 (^ Xt^~X^) coA. p A/^z^ -4 {5 ^ z<4 

•h SZ 0 |j^I’)j''X^ ^ <krR.p 6fl^Z.4 — CcfA© /<Am. 7 p ^ 
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5*Z. ^ ^ /4/<w 2. ^ ~h CfjrQ- 0 ^ C/T't 2^^ ^ "3 ® 

52 ^ I->j '" T-| J ^ -4-2 v. 0 /<iX^1.<\) ~ C^r<2, 0 ,^<jxav ^ /4xt^ 2 c|? 

+ 2 0 P 6dt!- 2 c|> ^ •+ ( 4^ "' 4^5^ 

~ ^ " A^vx,'^ ^ A^0A2 v ^ X/W'. <1? + D (^(j>) - <:<y^S^ A^<ly) 


(20) 


The preceding equations apply to arbitrarily large angular 
deflections of a rotating blade with no offset. 


2.2 Effect of Hinge Offset 

When an offset, e, is present, one introduces an additional 
axis system X,Y,Z which now represents the fixed axis system. 

The origin of the previous x,y,z system now circles about the Z 
axis at the offset distance e, while maintaining its axes parallel 
to the XYZ axes. 



Fixed axes — X X Z 

X = “X 4> e Cjr^S^ 

X* ^ G. /L-Vvs.^ 

= Z. 
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The velocities of the two axes systems are related by 

X — X - e. > 0 ^-^ ^ 

X = 1^ + e <><1.4- t 

2=4 


( 21 ) 


The kinetic energy T of the rotor blade is given by. 


T = i J [is I ^ 


vn 


( 22 ) 


+ 1^ -t xe + e. 03 _ ^ ^ d 


m 


This can be regrouped into the form, 

i 


T 


I (”r*^ / 

-\j?C + -‘-S } dm 

Jo 

e 

+• 6 $ 7 _ V f — nc ^ cIko 


(23) 


. ^ r 1 

+ 1 e ST. cl 
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The first terin above represents the kinetic energy of the blade 
found previously and given by (11) . The third term can be dis- 
carded since it is a constant, and will not contribute in Lagrange's 
equations (17) . The second term represents the additional kinetic 
energy T^ due to the hinge offset. This can be further simplified 
by taking = 0 in these equations to give, 



(24) 


Introducing y from (3) into the above gives. 


t; 


^ & c<x^ ^ 4 * ^ 51 ^ 

+■ ^ 0 ^ 4^ •t 0 dro. 51. ^ ^ 


(25) 


where is the static unbalance as in (16) . 

The additional terms in Lagrange's equations (17) coming from 
this T^ are. 


e 

E<]ua. I 

0 



j- 

fg ita. 1 

Sl 'e. ; 




Eguft. : 

z 

5^ e 
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These terms are to be added to the left hand sides of (18) , (19) 

and (20) respectively, and represent the sole contribution of the 
hinge offset. It can be seen from (21) and (14) that no additional 
gravity force terms will result from the hinge offset. 


2.3 Summary of Equations of Motion 

The nonlinear equations of motion given by (18) , (19) , (20) , 

(26) for the rotating blade can be simplified somewhat by expanding 
the trigonometric terms to 3rd order, i.e., 


<o 

^ i. - 


(27) 


and also by noting that for these blades, 




(28) 


Then considering only terms to 3rd order, the nonlinear equations 
(18), (19), (20), (26) reduce to, 

0 £~qu.g. - 

+ Z0^It^(5 + Z04>T^<|5 - 1^4 

+ 25z4'[i^ ■> I^e] I^e] 

+ Ag (e-9s)= 0<t.) 

( 29 ) 
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3 Ec^ucl : 


p + 0 + zefi^ 

-e''i^P P% aji4> 64>) 

+ [l-ji -|p*~ eS ec(>) + 

+ ^<^8 ( ?-?0 = t ® ■ r - ® t "- e T ) 

+ dS.-a.'P 



f r, - er^p - zep - e\+ 

-zj2.0[i^pcj> + r|e] ~ zs2.^i^(p+e<jj) 
+ sz'"[i„(gp+ eV + e$ 54 >] 

■e C^f {- 4> 4> } 


(31) 


Equations (29) , (30) , (31) represent the final nonlinear 3rd order 
equations to be investigated. It should be noted that in reducing 
the offset terras (26) to 3rd order quantities, the offset e itself 
v/as considered a first order quantity. 


I 
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SECTION 3 


LINEAR ANALYSIS OF FLAPPING-LAGGING ROTOR 

To gain understanding, the simpler flapping-lagging case will 
be examined before going to the full three degree of freedom case. 
A linear, small deflection analysis about a large static deflected 
position will first be conducted before examining the complete 
nonlinear analysis in the next section. 

The basic equations of motion for this flapping-lagging case 
are obtained from (30) and (31) by setting 0 and 0 equal to zero. 
The 0 itself is retained as a constant initial hinge setting for 
the blade. It will also be convenient to nondimensionalzie the 
time variable in the standard v;ay by introducing ip as, 

t ^ " (') = 

Under these assumptions the basic nondimensionalized equations 
become 

3 / 

OO o o o 

p(i- + ) - t- Z‘l’(^+e<p) 

* § ( 1 * e ■i-Vp*') + 0t = 

= + G, ( 9 - I-' - 8 I" - 0 i" ) 

E‘t|UgL > 

0 0 o , . 

<(> 4- 


(32) 


(33) 


(34) 
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where the follov;ing nondimens ional parameters have been introduced. 


e 


e i 


e«v^i/2_ 




mXV3 


— 

1 

n 7 - 

U) p 




'J€ 



1 


SI 


2 

- 



r>, SL"- \ ^ J 


3^ e 
^ £ 



UY\rfc>/tix\ 
lac^ e. y 


(35) 


2 .x 


i'oy' 


blp-elG. ^ 


In the G above, represents physically the natural frequency 

of the blade hanging as a pendulum v;ith no stiffness present. 

Equations (33) and (34) represent the basic nonlinear equations 
to be investigated for the flapping-lagging case. 


3.1 Static Solution 

A static solution to (33) and (34) can be obtained by neglecting 
all derivatives terms and by setting G = 0. Under these conditions, 
the equations reduce to. 


(i + e + v^-e''-|p-4>"')p + S't' = 


ep 


( 


£ -■ ■v'> + e*"- p'") 4> = 'j’i 


(36) 
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For a given S/ 0 , , B , v , v , these nonlinear equations can 

s s p (p 

be solved to give the static solution 4> = <()^, 3=6^* The 

solution is most easily accomplished by iteration, i.e., first 
2 2 

setting B and (J) equal zero in the parenthesis terms, then 

solving these linear simultaneous equations for B and (j>, then 

2 2 

correcting the parenthesis terms with the previous B and <{> , 
then solving for B and 4) again, etc. 


3.2 Small Perturbation Equations 

Having found the above static solutions B^ and (j)^, one can 
then investigate small perturbations of the nonlinear equations 
about these static positions by assuming solutions in the form, 

• <j> = ? 


(37) 


Here, B and 4> represent small perturbations about possibly large 
static positions B and 4> . Placing (37) into (33) and (34), 

O O Av rv 

retaining only linear terms in the perturbations B/ <j>/ arid cancelling 
out the previous static solution (36) , results in the linear 
equations , 




t- +• (i e + v-g -e -2-?!'- 1 ) P 

(e - ) ? - 

= [qe(i- - (6,ep„)p - 
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(38) 



o 


4* E c| u CL ; 
c> o 

? -z(g, + 8 4>„)p + 

> T O - ^„ ) <{> = 


/ _ 7. z- 

( e -H -f e 


1 


(39) 


These equations in B and (j) represent linear coupled equations 
with both forced excitation and parametric excitation present. 


3.3 Simple Linear Solutions 

One can obtain some simple linear solutions of the perturbation 

I 

equations (38) and (39) by arbitrarily uncoupling them. This will 
give an indication of the source and rough magnitude of B and (j) . 
Later a more accurate coupled linear analysis can be made. And 
of course, later still, in the next section, a complete nonlinear 
analysis will be made. 

A/ 

The B equation, v;hen uncoupled from the (p equation by setting 

/V 

all <p = 0 in (38) , is , 
oo 

" p.(i- ^ - ti) - (1 - r 1 
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Two types of strong oscillations are possible for this type of 
equation, a forced oscillation near ft and a parametric 

instability in the neighborhood of the first parametric instability 
region, ft % here, represents the natural frequency 

of this equation in flapping, namely. 


CO 


hi 


SI- 




(l e -4 


For the forced oscillation, one assumes approximately the 
steady state solution. 




-e t>, Qjro-'f' 


Placing into (40) and matching sin ip and cos \jj terms gives. 


a , 


(S(d (l 


2 . X 6 / 


1. 

o 


e -k- 


e + - s’- - z-p 


& 


The coefficients above can also be rewritten using the definitions 
of G and v_ in (35) as, 

p 


a. 


0 ^ z " 2 - c 

u)^ 






a ^ 6 

to 


4^: 




( 41 ) 


(42) 


(43) 


(44) 
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Equations (44) give simple estimates of the forced oscillation 
amplitudes in flapping. It is to be recalled that total amplitude. 


Ao= /~2 ! .2 . It can be seen from (44) , that resonance will 
P ^1 “l 

occur when the rotation speed reaches. 



— — ~ 1 

1 

N 0% 2.^^- e 


( 45 ) 


Since the denominator in (45) must be positive this resonance in 
flapping can only occur if 

0^ -h > e 


(46) 


This may v/ell occur for small hinge offsets e. 

For parametric instability in the first (and strongest) 
instability region, one assumes approximately the steady state 
solution, 

Q Ol, Ocn_ — (47) 

* 2 . « 2 - 


Placing into (40) and matching sin ^/2 and cos i|;/2 terms gives 
the equations. 


-o’" 

j^-“ G|9 b, = 



(48) 


O' I + 



■t 1 + e + 


0 




0 
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These homogeneous equations in and have the solution given 
by (47) only if their determinant equals zero, i.e., if, 






+- e 





( 49 ) 


0 


This equation can be solved for vt by bringing the second term 

p 

to the right hand side, then taking the square root of both 
sides, and rearranging to-give. 


2 . 

V* = 


-| _e 1-6% i (ep") 


(50) 


This gives two solutions corresponding to the tv/o boundaries of 
the first instability region. The boundary with + G is the one first 
reached here. Noting from (35) that G can be expressed as. 


(3 





(51) 


one may rearrange the criterion (50) to give, 




i io 


- _ z 


L 2- v'Z. 


M 




e" - i t 


+ 


T + (eP») 


-1 


(52) 


This defines the rotation speeds at which parametric resonance in 

2 

flapping occur. Since generally ^ windmills. 
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this can only occur if, 

0^ + A cb^ + > 2. -I- e (52A) 

which is generally impossible to meet. Hence no parametric 
resonance is expected in flapping. 

Looking next at the lagging case, the (}> equation, when 
uncoupled from the 6 equation by setting all 3 = 0 in (39) is. 


o o 

t 


+ ( e. 0 - j 4> = 

+ - G (t -€-%")? ] 


(53) 


The natural frequency of this equation in lagging is, 

60 , 


'N 


j e + - p.’- ‘ 


(54) 


and strong forced oscillations occur near fi w^, while the first 
parametric instability region occurs near ^ 'Z 2. U>^ - 

For the forced oscillation, one assumes approximately the 
steady state solution, 

<P ~ + (55) 


Placing this into (53) and carrying through as previously for the 
flapping 3 case gives. 
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L - ^ - 



Equation (56) give simple estimates of the forced oscillation 
amplitudes in lagging. It is recalled that the total amplitude, 


-r «2_ 


From (56) , resonance will occur when 


O4. 


Nj i- e + ^^-0' 


This resonance in lagging will occur if 

1 pj* > e + 0^ 

This condition is alv/ays met in practice, hence large lagging 
forced oscillations are distinctly possible at the rotation 
speed given by (57) . 

For pareimetric instability in the first instability region, 
one assumes approximately the steady state solution. 




H- 




% 

2- 


Placing into (53). and matching. sin ^/2 and cos ^/2 terms gives the 
equations. 


( 56 ) 


(57) 


(58) 


(59) 
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^ e + + s’-- - 1(1- - %")] . 

* [-|('^.-ep»)]^x =" ® ( 60 ) 

Ij i('^’o- 6pc.)J iij.+ [-^+a+‘'^4>+S"?o 


Setting the determinant of these homogeneous equations equal to 

n. 

zero as previously for the flapping 3 case gives. 


1 . 

= 


j_ 


2. 

— £ — 0 


Pc 



i. 2. 


(i- y'- % ) 


(61) 


Using the stability boundary ' above and expressing G as in (51) 
previously, (but with reference to co^ now), the criterion (61) 
may be rearranged to give the rotation speeds for parametric 
resonance as, 





( 6 . 2 ) 


. 2 

Since generally ^ windmills, 

resonance will occur if. 


this parametric 


_l_ 

4 


+ 


X. 




> e + 


(63) 
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This condition is generally met for small offsets e, hence 
parametric resonance in lagging is distinctly possible at the 
rotation speeds given by (62). 


To end this brief discussion of simple linear solutions, it 

may be of interest to introduce some simple physical interpretations 

to some of the results, equations (41) to (63). The criteria for 

strong forced oscillations in flapping (46) and lagging (58) , and 

the rotation speeds at which these forced oscillations occur (45) 

and (57) , come simply from the requirement This can 

easily be shown by setting n = in (41) and (54) respectively. 

Similarly, the criterion for the first parametric instability 

region in flapping (52A) and lagging (63) can be shown to come 

from the requirement = 20)^^. The rotation speeds at which these 

parametric instability regions first occur (52) and (62) are 

2 2 

slightly less because of the (Wp^ ^‘^Pend^^cJ)^ factors. 

If ^ ^ ^ ^‘^Pend'^*^4>^ then these rotation speeds 

also occur at = 2o)„. 

N 

Another quantity can be given a simple physical interpretation, 

2 

namely, the ^“pend'^^g^ factor appearing in (44) and (52) . Consider 
a blade held horizontally against gravity and restrained at the 
3 hinge by the spring as shown belov;. 

p 



The angular deflection 6^^ under its own dead v/eight is given by. 
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( 64 ) 



P 


where dm is the static unbalance about the n axis. The 

natural frequency Wq of this blade is defined as 

p 


CO 





(65) 


Solving (65) for k., placing it into (64), using the definition 
2 ^ 

(w„ j) =g S^/I given in (35), and noting that S , results 

Pend ^ ? n ^ 5 q 

finally in 


e = 

1 

(66) 

V DVY 

V Wp / 

2 

Thus, the factor 

simply represents the 

static deflection. 

B of a horizontal blade under 

its ov/n dead weight. 

Similarly it 

can be shov/n that, 

^ =. 

.7- 

(67) 

•jpw 

V t0<|, / 


These and give simple interpretations of the forced 
amplitude oscillations (44) and (56) . 

In summary, the simple linear analysis in this section has 
indicated that gravity effects are likely to be more important 
for lagging motions than for flapping motions, and that both strong 
forced oscillations and parametric instabilities are possible for 
lagging at the rotation speeds. given by (57) and (62). 
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3.4 Complete Linear Solution 


complete linear solution of the coupled small perturbation 
equations (38) and (39) can be made for more accuracy. For the 
case of forced oscillations, one would assume the approximate 
steady state solution. 


AX 

? 







j ^ +• la j C-Cf^ ^ 

Cl-2, 


( 68 ) 


The higher harmonics, sin2t/j, cos2i(; will be neglected here. 

Then placing these into (38) and (39) and matching the constant, 

the sin ip, and the cos ip terms of each equation results in six 

equations in six unknowns, $ , a,, b. , (J> , a», b_. The 0 and 

O X J. C ^ o 

<() are included in (6S) to allow for small center shifts from the 

static values 0^ and The solution (62) will nov; include the 

effects of the small second parametric instability region near 

^ 0 )., in addition to the dominant forced oscillation resonance 
N 

there. The six linear equations can be readily solved by inversion 
to give the forced oscillation amplitudes in (63). Also, the 
determinant of these equations can be evaluated numerically for 
different values of rotation speed Q to find what makes the 
determinant equal zero. This would then represent the boundary 
of the small second instability region. 


Instead of developing and presenting these coupled linear 
equations here, it will be more convenient to present and solve 
them later as a subcase of the complete nonlinear equations to 
be given in the next section. 

For the case of parametric instability in the first instability 
region, which is always the strongest instability region, (see 
Bolotin [40]}, one would assume the approximate solution. 
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(68A) 




di 




■h 


b, Ccrfi- — 


cp 


- 

a, i +- 

^ 1. 


k 

i i, 


and obtain the appropriate four equations in the four unknowns. 
The higher harmonics sin3i|>/2, cos3t|^/2 are neglected here. The 
determinant would then be examined numerically to find what 
rotation speeds Q make it equal zero, thereby determining the 
stability boundaries. Again, this linear solution will be 
obtained as a subcase of the complete nonlinear equations in 
the next section. 
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SECTION 4 

NONLINEAR ANALYSIS OF FLAPPING-LAGGING ROTOR 

The linear solutions given in the previous section serve 
only as a guide to the small amplitude behavior of the rotor 
blade. When the amplitudes become large, the nonlinear terms 
in the equations will serve to limit the predictions of linear 
theory. This is particularly true of the forced oscillation 
resonances and the parametric instability regions where infinite 
amplitudes are predicted by linear theory. Accordingly, the 
complete nonlinear equations (33) and (34) will nov; be examined. 


4.1 Forced Oscillations 

For forced oscillations of equations (33) and (34) , one seeks 
limit cycle solutions in the form 


^ ^ 4- a, -4- b, 

The higher harmonics, sin2if;, cos2(|/, etc. v/ill be neglected here. 
The 3^ and ^ above represent the total centershift from zero 
and would now include the static solutions 3^ and c})^ plus any 
additional centershift due to the oscillation amplitudes, a^ bj 
b^. Placing these expressions (69) into (33) and (34) and 
matching the constant, sin \p , and cos \jj terras of each equation 
and discarding the higher harmonic terms, gives after much algebra 
and trigonometric reduction, the following six nonlinear equations, 


(69) 
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F« 



- 


- 



K 

1. J 




= 0 


( 70 ) 
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where the elements F. . 

i3 


are given as, 
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+ e - <C- 
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^31= Q - NG, I j(a?-+sb:) tK + sb,")] 

f"3:u= 

£ - g - - N 5 i (aN k^; - t 

Fj^ = - N Q j :^ ( a.a^ + 3 U, b».) j 

F35- = 2, ( 0 4^^) 

Fjt = 0 - zp.t. 1 

F j = 0 *" N ^ (Z,Q.jb b, b>3, ^ 

F^^ = G,i6t 

F^3 = -G, - NGj [ -^a,a^^ 

F^^ = 2 ^ e + e" - p," 

= -^i(t.-epc) 

F^, = G,(i-g--fc")- NlG,[fJaN3br)-t(«t+br)j 
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NGi 

0 - - N I b, bi."^ 

z(?. + e+^) 

-Nc^ 

2j + e- !■*■ 0 bJ^-^(a 

0 

= - N Ci ^ t + 3 k, k)} 

^2. = - 

Ffc5 = 9 - - N 

f^4-= g, (i-4"-^) " 

F,r - 0 

-F(,i, = 2, + e - b + e"" - ?c + N I i ( « 
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Fs3 = 

F5. = ■ 
Fs. = 
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b,^;] 


(a^ 3 bt)j- 
^- bf)} 
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f;o 

= - £ fj. 

- N 

ai-b, ] 

+ g(3a):+b3]| 


= -qe(i- 

s!: - p" _ 't’-N + 

(» t: t J ^ 

N Q ^ 0 g (,? bO ' 

^30 

= NG, ^ e 

^ (a,b, -t a 2 .kj. 

)] 



= - 

N ^ 0 (a,bj_ 

- <l2.b,)] 



= 


- N(q 


F,o 

= - N (q ^ 

- ie(» 

,,bj + 



The nonlinear 

terms appearing 

in the F. . expressions 

have been 


grouped to provide a reasonable symmetry. Other groupings could 
also have been chosen. 

In the above elements the following definitions have been 
introduced for convenience, 



( 71 ) 
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For a given configuration defined by the six parameters 9,3/ 

— — 2 ^ 
(j)^, e, ojg, and “ *^DW' solves the nonlinear 

equations (70) for various values of rotation frequency ratio 

fi/o).. In evaluating the coefficients, it is convenient to use 

the relations. 



The solution of the nonlinear equations (70) is best accomplished 
numerically by using an iterative Newton- Raphson technique, v/hich 
uses some initial estimate of the solution to begin the process. 

As a start one may use the trial solution B = ^ , 4>^ = <l> r 
a^ = b^ = a^ = b 2 = 0. One may be further guided by the simple 
linear solutions (56) and (44) given in Section 3.3 However, one 
must realize that for nonlinear equations, multiple solutions 
may be found in certain cases. The particular solution obtained 
would depend on the initial estimate in the iterative computation 
process. Solutions would probably vary W’ith frequency as shovm 
below. 



The effects of various parameters such as 3^/ 9 etc. could 

be assessed as desired. Also, the complete linear case mentioned 
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in the previous section 3-h could be worked out by simply setting 
N = 0. 

Finally, it should be mentioned that some of the limit cycle 
solutions obtained v;ill be unstable and as such, have no physical 
reality. These unstable solutions can be checked formally by a 
stability analysis or may be inferred by experience. 


4.2 Parametric Resonance 

For parametric instability in the first instability region, 
one assumes a limit cycle solution of (33) and (34) in the form, 


4* 4'o ■*" +- 

The higher harmonics sini|j, sin3i(;/2, etc. will be neglected here. 
Placing these into (33) and (34) and matching the constant, sin i|i/2 
and cos ^/2 terms of each equation and discarding the higher harmon 
gives again six nonlinear equations. 









r ^ 

r ^ 

R, 

Pii 

P,a 

P 

r j4_ 


F.. 


Pc 

1 

p 

*10 

P., 

p« 

P.a 

* 

T 

■ - 


a, 



Pa. 

Pa. 



* 

• 

> 




P4-1 




* 

9 





P^, 

•9 


% 

- 

9 


Oz. 



1 


- 


t 

P 






(73) 


ics , 


( 74 ) 


44 


where the elements P^j are given as. 
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In these elements, the same definitions (71) and (72) are used 
as before. The equations (74) are solved as previously for the 
forced oscillation case, only now, one is guided by the simple 
linear parametric instability solutions (52) and (62), and one 
seeks solutions near 2. These steady state limit cycles, 

appearing at roughly twice the rotational speeds R of the forced 
oscillation resonances, may be equally as severe as the forced 
oscillations. It is to be noted that the actual vibrations them- 
selves occur at roughly w ^ aK even though the rotational speed 

(p 

here is roughly R 2w,. 

<{> 

Other nonlinear subharmonic and superharmonic solutions can 
be investigated in a similar manner by introducing equations (69) 
or (73) v/ith additional harmonic terms present ii.to the basic non- 
linear equations (33) and (34). One would then harmonically 
balance these additional harmonics v/hich would lead to larger 
size nonlinear algebraic equations in place of (70) and (74) . 

Such subharmonic and superharmonic solutions for simple beams 
were examined by Tseng and Dugundji [38]. 
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SECTION 5 

NUMERICAL RESULTS FOR NONLINEAR FLAPPING-LAGGING ROTOR 


Numerical results using the previous nonlinear analysis of 
Section 3 are obtained for the following configurations: 


CASE 
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. 088 
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flexible 

lag rotor, i 


the case for which oj,= 3.37 ca , 

(j) pend . 

5.1 Forced Oscillations 

First, simple uncoupled solutions are calculated for the 
different cases using Eqs . (36) and (44). Theuthe linear and 

nonlinear solutions are obtained from Eq. (70) by computer. For 
the corapution of nonlinear solutions, the corresponding linear 
solutions are used as initial guess as far as possible. The 
results for the above three cases arc presented in Tables 1-3. 

From these results, one can see that the simple solutions give 
reasonable estimate of linear solutions, particularly, where 
coupling motion is weak. Also, linear solutions agree v;ell with 
nonlinear solutions away from the resonance region. This is 
quite apparent from the fact that the linear theory is good for 
small amplitudes and becomes inadequate for large amplitudes 
which take place near resonance condition. All these solutions 
are checked for stability by giving linear perturbations to the 
steady solution and then studying the growth of these perturbations 
(discussed later in Section 5.3} 

In Fig. 1, the solutions are plotted for Case I. For this 
configuration of zero coning angle, the response is uncoupled and 
one gets only lagging amplitude. Here it can be seen that the 
linear and nonlinear solutions are quite close except near resonance 
(i.e., ~ .95 to 1.15). Arrow mark on the graph is a resonance 
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point obtained by linear theory. The nonlinear solutions become 
unstable for A, greater than about 1. 

In Figs. 2a, 2b and 2c are respectively plotted lagging 
amplitude, flapping amplitude and center shift for coupled 
response of Case II. Here, one can find a large shift of non- 
linear behavior from the linear one particularly after the ampli- 
tude started increasing. At large amplitudes, the nonlinear 
resonance curves bend tov/ards decreasing frequencies depicting 
a typical nonlinear softening spring type system. From Fig. 2a, 
one finds that for lagging amplitude, the higher frequency branch 
becomes almost flat in the overhang whereas lower frequency branch 
becomes unstable after maximum slope. Similar strong nonlinear 
characteristics are also visible for flapping motion in Fig. 2b. 
Further it can be seen in Fig. 2c that the mean angle setting 
about which limit cycle oscillations take place increases with 
increasing amplitude. These curves show that small initial coning 
angle of the order of 9° can produce an appreciable change in 
the nonlinear forced response of the blade. 

In Figs. 3a to 3c, the results are plotted for Case III. 

This physically signifies a system in which lagging hinge stiff- 
ness is lov/er than that of flapping. Here, also, like Case II, 
one sees the softening spring characteristics for large amplitudes. 
In fact, the nonlinear response for case III appear to be more 
violent than Case II. 

5.2 Parametric Resonance 

The numerical results are obtained by solving Eq. (74) by 
computer, for the same configurations for which forced oscillation 
response results were obtained. These results are plotted in 
Figs. 4-6. 

In Fig. 4, corresponding to Case I, results are plotted 
for linear and nonlinear solutions. Like forced response of 
Case I, (i.e. zero coning angle) one gets here also only uncoupled 
lagging motion. It is found that there is distinct unstable band 
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for of 2.52 to 2.64. The linear solution predicts infinite 

<P 

amplitude in this band whereas nonlinear solutions give two 
limited amplitude branches which try to bend toward each other 
in this region. However, for higher lagging amplitude (-1.4), 
the nonlinear solutions become unstable. The trivial solution 
of 3(,“^2.~^l~'^c~^2~^2~*^ exists for every except for this 

instability band where trivial solutions becomes unstable. 

In Figs. 5a to 5c, results are plotted for Case II. With 
a small coning angle (-9°), there is hardly any change in the 
linear response but there is distinct change in the nonlinear 
parametric response of the blade. Nonlinear response is coupled 
one with lagging motion more dominant as compared to flapping 
motion. Like forced oscillations of Case II, the nonlinear 
parametric resonance curves for large amplitudes bend towards 
decreasing frequencies, again, depicting a typical "Softening 
Springs" type system. The higher frequency branch for lagging 
motion, becomes flat and extends to lower frequency region 
and lower frequency branch becomes completely unstable. Again 
the trivial solution is stable except for of 2.52 to 2.64. 

In Figs. 6a to 6c, the results for Case III are plotted. 

The results are just similar to those plotted for Case in II in 
Figs. 5. Again one gets the coupled response with strong lagging 
motion. 

5.3 Stability Analysis 

The nonlinear solutions obtained by Harmonic Balance method 
for forced oscillations as well as parametric excitations are 
not always physically existent. One has to make stability check 
on these solutions to find out whether any of these can be a 
physical reality. So, these solutions are further investigated 
here for stability by giving small perturbations to the steady 
solutions and studying the growth rate of these disturbances under 
the assumption of slov;ly changing functions. The solution will 
be unstable if the growth rate of perturbations with time is 
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positive. By slowly changing is meant here that the increase 
of a function during a period is small as compared to the average 
value for this period, i.e., 

o 

^ 2FT 

a. 

See Bolotin [40] for details. 

First, considering the stability of forced oscillations, 
one can write perturbed solution in the form 




<< 1 


do 

0 


zjr 


<< 1 etc 


(75) 


where 6-^^, ^lO'^lO' '^cO'^20'^20 ^®P^®^ent the steady-state 
solution for which stability check is being made, and 3 ,a,,b , 

/V ^ ^ C J. X 

*^c'^2'^2 time dependent perturbations given to respective 

steady components. To study the growth of these disturbances, 
the above equations are substituted in the basic flapping and 
lagging equations (Eqs. 33 and 34). Then retaining only linear 
terms in perturbations and their first order derivatives and 
subtracting the steady-state solution ,once again on matching the 
constant, sinif/ and cosip terms from these tv/o equations, one gets 
six linear algebraic equations which can be concisely put into 
matrix form as 
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( 77 ) 


Putting , 


Substituting this into the above Eq. (77) , results in an 
algebraic eigenvalue problem which can be easily solved using 
any standard eigenvalue subroutine. The solution will give six 
eigenvalues — if any of the eigenvalues has positive real part, 
this means perturbation v;ill grow with time and this marks the 
solution as an unstable one. 

On similar lines, the stability of parametric excitation 
solutions are checked by putting, 

P = ^ r«,o+ fi+S('.^)Jc»sar 

'P ~ [fc + t(ip)]cosv 


Following the same steps like Forcing Solution v/e get here also 
an algebraic eigenvalue matrices. 



Again the nature of the roots X will indicate whether the solution 
is stable or not. 
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PART B: EFFECT OF AERODYNAfilC FORCES 

SECTION 6 

FORMULATION OF AERODYNAMIC FORCES 


The aerodynamic forces are obtained using quasi-steady 
airfoil theory. The elemental lift and drag forces acting on 
airfoil , section can be written. 


tlL = if\V a c c(»>. 

Co 

v;here 

a = section lift curve slope, 

= blade profite drag coefficient 
c = blade chord 
p = air density 

= resultant velocity, /y2^^2 

6^. = local blade built-in incidence 
= inflow angle, tan ^ ^ 



( 81 ) 


[r is same as Cl 



As shown in the above figure, U,v are flow velocity components 
along blade axes n and ^ and 0^ -a^ is the effective angle of incidenc 
Resolving the aerodynamic forces along the blade axes, one gets 


dN = dL c^tTt - cfD Mac5d^ 
dC = dL Mnoc. -H dj) co'iPC^ 


( 82 ) 



Assuming the effective angle of incidence is small (i.e., 
below stalling angle ) , such that 

$l~ (Ki 


and also considering j 
written as 


1 

^ u2 


<<1, the forces dN and dC can be 


clN « f a c cLj^ [ U^sia®. - v Uco^>®^ - ^ \j l/J 

dC 21 f a c CcU V L/ 5ia 0^ - V^( cos (9^ - a ^)] 


(83) 


From Eq. (10) of part A, the angular blade velocities about the 
three axes r\/ ? are given as, . 

6J^ = ^(skS s(a4’ +• c.os<9 slap cob<^)~ ^ siac^ 

=. <n.( sia (9 cos4> - cosQ siap Sincjp) - (I co!?c|> (84) 



The inflow velocity through the rotor is, 
U- ^ -'Kri a 

v/hich can be resolved along three blade axes as 

Hi = + Ut'2 tj. 


(85) 


(85a) 
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where 


= -ylARCcosG siap cos + sia 0 sia cjj) 

= -\sik(~coiQ siap sia4 +- sia0 coscji) 

-X^R COS0 oosp 

In the above, X is the inflow ratio parameter and is assumed 

here constant over the disk area for a uniform approaching wind stream.* 
Introducing the blade hinge-offset effect, the velocity of 
the hinge-point is given as 

S, (86) 

which can again be resolved along the three blade axes as, 

\ ^ ^ 

where 

- en.(- staG siap cos^ + cos0 s(a<4) 

smG sva^ <^s6 coscj?") 

= _ e^l Sia 6 005 p. 

The resultant flow velocities U and v for any point 5 on the 
blade is then, 


XJ ~ h. CO^ — 


4- LL, 


'I 


- >i(^Cos0 cosp + c^) + sin<4 +s!n0 cos^) 

+ Ciil(SuriS Sia"^ S>r\<4 + cos © Cos <4^ 


(87) 


V ^ ^ Uic^ + 


= COS0 

~ CXI 5ia 0 


S\tip> sm4 - sia© coscb) + ?iS cos^ + AaR ca5Sco>^ 

( 88 ) 


go's p. 


This is strict] y true only for an ideally twisted rotor. For 
other 1 ‘otors, there may be some variation with r. 
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Like part A (Eq. 27), expanding the trignometric terms to third 
order in above equations we get 


U = ah.(,- 4 


(89) 


with the inclusion of aerodynamic forces, the generalized 
forces Q_ and Q, in Lagrange's equations (Eq. 17) have the 

p ? 

additional contributions, 

■ (90) 

sjjan. 


One may assume the built-in incidence along the blade span is 

where 0^^ is constant incidence along blade length and alone 
represents the ideal twist distribution. 

After performing various integrations and combining with 
Eqs. 33 and 34 we get the nonlinear equations of motion for a 
flapping-lagging rigid rotor with f lexible^springs at the hinge- 
point in the presence of airflow as 
Flapping Equation; 


? 0 - 4 ) - 1+ e + 

-^4 ^ Cp Vp P - I 4 G(-® + f' 4 4 e£)slnr^ 


= ^ 4 4- c,^4^ 


P 4 C, 2 . + 4 ^15 ^44 ^4 -+ c,s -h C6, 4 cp^qj 


(92) 
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Lagging Equation: 


$ + I + t ( ^ + >4) + <5^ H- Q^4) - + g c| 

-VcJts + G('-f +Q^4>)smH' + ^‘)cosy 

= -X [B^ ^ B5 + S, t ^ 4- !> aj -h p 8 ,^ + p fB„ 

o 

"■ P ^ ^3 t- 6,4-^ S,s 4-p<^q + 4>^^o _ f BC,J 

( 93 ) 


In above equations, y = pacR /I^ is the Lock number and, Cg ^nd 
represent structural damping coefficients for flapping and 


lagging motions, respectively, 
in Appendix I. 


The other constants are defined 



SECTION 7 

NONLINEAR RESPONSE OF FLAPPING-LAGGING ROTOR 
WITH AERODYNAMIC FORCES 

Seeking a general limit cycle solution of nonlinear Eqs. 92 
and 93 for both forced oscillations as well as parametric resonance 
in the form of 

■p « I + a, sia pcH' -v \ cos ocH' 

-t aj^sinccv 4- \>^cosoc^ 

where a is the ratio of response frequency to forcing frequency. 
Here a=l represents forced oscillations and ct = gives para- 
metric resonance. Substituting these expressions (94) into 
Eqs. (92) and (93) and balancing out the constant, sin a\p and 
cos aip terms of each equation and neglecting higher harmonics, 
gives six nonlinear algebraic equations. 



The various terms in the above matrices are defined in Appendix II. 

The solution of these nonlinear algebraic equations (95) can 
be again obtained numerically by Newton-Raphson technique. 
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stability Check of Solution 

As discussed in Part A, the nonlinear solution got by 
Harmonic Balance method has to be checked for stability to prove 
its well-posedness . Once again, here the solutions are checked 
for stability by giving small perturbations to these steady 
solutions and studying the growth rate of these perturbations 
with time under the assumptions of slov/ly changing functions. 

(See Section 5.3 for more details.) 

Writing the perturbed solution as 

^ + [Q(o a,(‘t)]sin.cciK 4 - [b>^ + t,(v)] COS0C4/ 

(96) 


where 3 n, a,-, b,,,, <p a^,., b-,. is steady solution for which 
cO 10 10 ^cO 20 20 ^ ^ /V ^ 

Stability check is being made and 3^, a^, b^, are the 

time dependent perturbations given to respective steady solution 
components. 

These equations (96) are put into basic nonlinear equations 
of motion (92) and (93) . Retaining only the linear terms in 
perturbations and their first order derivatives and filtering 
out the steady solution, then again on matching the constant, 
sin and cos a\p, one gets six linear algebraic equations. 
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writing 
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( 98 ) 


substituting (98) into (97) results into an algebraic eigenvalue 
problem which can be solved by using one of the standard subroutines. 
Solution will give six eigenvalues X and the nature of the eigen- 
values will decide whether solution is stable or not. If any of 
the these eigenvalues has positive real part means perturbation 
will grow v/ith time and makes the solution unstable. 


Numerical Results 

Numerical calculations for forced oscillations as well as 
parametric resonance with aerodynamic forces are repeated for 
most of the configurations of Part A. Once again 
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For aerodynamic characteristics. Lock number y of 12 is 

used for most of the results and the values of 'a' and are 

d 

taken respectively of 6.0 and .012. 

7.1 Forced Oscillations Resvilts (a=l) 

In Fig. 7 results for Case I of zero initial coning 
angle with aerodynamic forces are plotted. Here X=0 represents 


\ 
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a configuration with zero inflov; velocity. It can be seen that 
for this case of no inflov;> response is uncoupled and. one gets 
only lagging amplitude, which is quite similar to that of no 
aerodynamic case (Fig. 1). Here also the amplitude shoots up 
near resonance (— - 1.05). However, v;ith the inclusion of small 
inflow ratio of X-.l, the amplitude is well bounded. Also, with 
inflow, response is coupled, though the amplitude of flapping 
motion is comparatively much smaller except near resonance (Fig. 7b) . 
Further, in Fig. 7b, the mean angle about which limit cycle 
oscillations take place, is plotted and has negative value 
because of thrust direction for this positive inflow. 

In Fig. 8, the lagging amplitude response for Case II 
with aerodynamic forces is plotted. Here positive X represents 
configuration with initial coning angle facing into flow direction 
and negative X the case of initial coning angle facing along the 
flow direction. For X=0, the response amplitude is entirely 
different from that of no aerodynamic case. Fig. 2a. With the 
presence of aerodynamic forces, the instability overhang region 
completely vanishes. Though the amplitude increases near — ^ 1.05, 
but it is bounded with flat peak. The peak value is higher^for 
negative X as compared to positive X. This can be explained with 
the help of Fig. 7b that with positive X, the coupling is reduced 
by decreased B and for negative X the coupling is amplified by 
increased B^. It is found that the flapping amplitude response 
though smaller in comparison to lagging amplitude, the peak for 
negative X is quite higher than that for positive X (not shown in 
figure). Further, it is seen that the solution for X = -0:1 
becomes unstable in the peak value region. 

In Fig. 9, the results for Case III are plotted. For 
this configuration of softer lagging hinge stiffness than flapping 
stiffness, the behavior is quite similar to that of Case II in 
eliminating the violent overhang instability region with the 
presence of aerodynamic forces. Here the peaks for lagging 
amplitude are comparatively sharp but are of nearly same magnitude 
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for the three different cases of X=-.l, 0, .1. This is mainly- 
because of higher flapping stiffness, the average coning angle 3 
is not too much effected by the small value of inflow, X=+.l. 

Figure 10 illustrates the lagging response for a typical 
configuration of Case II with X=-0.1 for different Lock numbers, 

Y- Here, one can see that by reducing Lock nximber, the behavior 
of forcing response tends towards that of no aerodynamic case 
(Fig. 2) . For Lock number of 3, there is little overhang tov;ards 
left depicting a typical nonlinear softening effect. Again, the 
higher frequency branch remains stable in overhang region whereas 
lower frequency branch becomes unstable. 

In Fig. 11, the effect of structural damping on forced 
oscillations with no aerodynamic forces for Case II is presented. 

For small value of damping coefficients ?g=C^=.01 encloses the 
stability branches (see for comparison Fig. 2a) . The behavior 
is quite similar to that shown in Fig.. 10 for low Lock number. 

This means that the aerodynamic forces can be looked as if they 
add the equivalent structural damping to the system depending 
upon Lock number y inflov/ X. 

7.2 Parametric Resonance Results (a=-|-) 

In Fig. 12, the results for Parametric Resonance of 

Case I in the presence of aerodynamic forces are plotted. There 

is a distinct band of instability for zero inflow and the response 

is uncoupled lagging motion. From stability analysis one can see 

that the trivial solution of 3^=a,=b,=c^ =a_=b_=0 is stable for 

C i J- c ^ 

every’— except in the instability band where trivial solution 
becoraes'^uns table. This behavior is quite similar to that for no 
aerodynamic forces (Fig. 4) . The two instability branches bend 
towards each other and the higher frequency branch also becomes 
unstable for large lagging amplitude (A^~.75). However, with 
the inclusion of small inflow X=.05, this instability band 
completely vanishes. 

Figure 13 presents the response for Case II. With the 
inclusion of aerodynamic forces with no inflov/, one can see that 
the behavior of rotor is quite different from that of no aerodynamic 
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solution (Fig. 5) . The overhang instability region gets comple-tely 
eliminated and the tv/o stable solution branches joins in the in- 
stability band resulting in the limit cycle amplitude. The 
response is coupled but the flapping amplitude is of an order of 
magnitude lower than the lagging amplitude. Again with the inclusion 
of small inflow, the instability band vanishes. 

Figure 14 shows the effect of Lock number on Parametric response 
of Case II with A=-0.1. For Lock number y= 12, there is no instability 
region, but for and 3 there is instability region with overhang 

towards left examplifying a typical nonlinear softening effect. Out 
of the two branches of solution, the higher frequency branch is stable 
while the other branch is unstable from the very initiation. The 
other noticeable thing about 'these graphs is that with the increase 
of Lock number, not only the instability overhang and peak amplitude 
reduces but also the instability band shifts towards the right. 

Figure 15 illustrates the effect of structural damping on 
parametric oscillations of Case II with no aerodynamic forces. 

Here the solution is not very different from that of no-damping 
case, Fig. 5a, except that the two branches of solutions shrinks 
in the overhang region. 


SECTION 8 

SELF-EXCITED FLUTTER RESPONSE OF FLAPPING- LAGGING ROTOR 
In the earlier sections it has been shown that a rigid rotor 
restrained by two springs at hinge point representing flapping 
and lagging motions can get into large nonlinear response due to 
gravity forcing excitation. Here, it is intended to discuss that 
this torsionally rigid rotor, in the absence of gravity forces 
can also lead to self-sustained oscillations caused by the inter- 
action of aerodynamic forces with structural vibrations of 
blade. 

The equations of motion for flapping-lagging rotor for flutter 

can be rewritten from Eqs. (92) and (93) by eliminating gravitational 
effects. 

Flapping : 

p(l-<(5^) - 2^ ^ ^ ^ e<p) 

^ ^ 1^0 "** ^7 ■I' P ^ ^10 

(99) 
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3.1 Simple Linear Solution 

It is sometimes advantageous to deal with simple solutions to 
get the feel of the nature of phenomena and then later on using 
these simple solutions as a guide one can make more rigorous 
analysis. Here also first simple linear solution for self-sustained 
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oscillation is worked out by assuming that the small amplitude ' 
motion is taking place about possible large static positions 
Bq and 4>o* 

P = I + ? 

c{> = c|) -V c|> (101) 

where 3 and ({i represents small perturbations. Placing (101) into 
Eqs. (99) and (100) and keeping only linear terms in perturbations 
and their derivatives and subtracting the static solution results 
in 


3 Equation: 

<50 

t('^8 +1 to '^10 + + 4 C,j)h-|(C, Qs)] 


(102) 


(j) Equation: 

aa oo 2 2.2 

= - I [ pC’Po Bs + 1 S,) + ? t S4 + I S 7 > ^ s t ^0 l"s„ + 4 B,,) 

+ $(B,+ It6,3+ " 

(103) 


and static solution equations are 

1 ( 1 + e A- )^-e -f-l-4) +t6='^*| + ^'fCoA.g\ A.<|q 

|s * t(^ * 4 ’' + 

^ (104) 

These equations are further simplified by considering a rotor 
with ideal twist distribution (only B and with no hinge offset 
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(e=0) and also assuming that the static equilibrium angles 3 q<1 

2 

and so that and (J)^ can be neglected. Equations 102- 

104 reduces to 


3 Equation; 


00 o ^ \ 2 


(105) 


(p Equation: 

00 ^ ^ 

? + + +5S^>')+5i;^y^]+(/-^p^)f (106) 


and static equilibriiam equation 

i) - ih * 

For the critical condition of flutter 



(107) 


(108) 


substituting in Eqs. (105) and (106), gives characteristic equation 


(tco) 4- too VYl. -V m 

p p 

ico m* + wIa 

P 


ICO n. 
P 

(iwt + ico + n^j, 


_ 


. j 


(109) 


For non-trivial solution, expanding the determinate and comparing 
the real and imaginary parts separately to zero, one gets 


Real part; 

Imaginary part; 

'■(y>^ + ) - (-n^ "'p ”.p - 


CO 


( 110 ) 


( 111 ) 
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where 


= 5P-x?t + xA. 

Y To 3 g 

n. = -5^ + 

P * 6 3 

I 4 . =i(lC^+2l9^A)j-^!'^4 
? " |'7^l[f(®2-'') + ^’■|] 

Equations 110 and 111 can be solved for flutter boundary by 
simple iterating scheme. 

8.2 Nonlinear Solution 

Assuming the limit cycle solution of nonlinear Equations (99) 
and (100) for self-sustained oscillations of the form 

^ ^ ^ 4 a, sinoc^f + b, cos 

c|> ~ -t a^sin ^cH/ cosocif 

where a represents the ratio of flutter frequency to rotational 
frequency. 

Placing these expressions (112) into Eqs. (99) and (100) and 
matching the constant, sin ai|; and cos a4> terms of each equations 
and neglecting the higher harmonics yields six nonlinear algebraic 
equations. Since v/e are interested in the limit cycle flutter 
solution, this can be done by assuming a known lagging amplitude 
of flutter and finding the other associated unknowns to get this 
condition. Thus putting t>2=0, and which now represents lagging 
amplitude, a known quantity, into these nonlinear algebraic equa- 
tions, 
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(113) 


2 

where Z = (w./fi) and the other terms in the above matrices are 

4 > 

defined in Appendix III. 

The solution of these nonlinear Equations (113) are obtained 
numerically by iterative Newton-Raphson technique. 


8.3 Stability Check of Solution 

The nonlinear solution got above using Harmonic Balance 
method is checked for stability again by giving small perturba- 
tions to the steady solution 


P = Since t,(.M-')] i^oScCH' 

^ " tco + + [Qjo + 


(114) 


v;here 3 ^q/ a^, b^, 4 >^q/ ^2 with corresponding Z and a 

is the steady solution and a^, 6^, ^ 2 , ^2 with same Z 

and cc are the time dependent perturbations. 

Following the same steps like forcing response, placing (114) 
into Eqs. (99) and (100), keeping only linear terms iii perturba- 
tions and subtracting the steady solution, once again on matching 
the constant, sin a'jj and cos arranging appropriately one gets 
in the form of algebraic eigenvalue 
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(115) 


This is just the same as Eq. (97) . Again, the nature of roots X 
will decide whether the solution is stable or not. If any one of 
the roots has positive real part, points out that the solution is 
unstable. 
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8.4 Numerical Results 

For most of the flutter calculations, Lock number of 12 is 
used and values of a and are taken respectively of 6.0 and 

.012. It is also assumed here for flutter calculations that 
initial lagging angle and initial built-in incidence are zero. 

X=0.1 is taken for most of the results as this inflow represents 
nearly the maximum power extraction condition of the rotor. The 
critical flutter boundary is obtained from nonlinear analysis for 
very small lagging amplitude of .05. 

Figure 16 presents the critical flutter boundary for centrally 
hinged blade v;ith no preconing (i.e. e=0. 3^=0) for different 
inflow X. Both, simple linear solution and more accurate Harmonic 
Balance solution are plotted for comparison. The region of 
instability lies inside the respective contours. It can be seen 
that simple solution gives reasonable estimate of flutter boundary 
for higher values of lagging stiffness, v^, and it starts deviat- 
ing more and more at low v^. This is, mainly because in simple 
solution, cJ)q, the equilibrium lagging angle is neglected and it 
becomes appreciable at low values of v^. One can also find that 
the instability envelop increases in size with increasing inflow 
ratio X. This can be explained by looking at linear Eqs. (105) 
and (106) that X effects the potential destabilizing flap-lag 
coupling terms. 

Figure 17 shows the effect of preconing and hinge offset on 
the flutter solution of rotor with a typical inflow X of 0.1. On 
comparison with the corresponding instability graph for X=0.1 from 
Fig. 16, one can easily visualize that with the hinge offset the 
instability envelop not only expands in size but also opens up at 
lower end from the elliptic shape. Preconing also has an important 
effect on the flutter boundary - positive preconing, coning facing 
flow direction shrinks instability envelop v^/hereas negative pre- 
coning expands the instability envelop. This is because the blade 
preconing effects blade equilibrium coning angle and thereby effect 
mainly the cross coupling terms. 


I 
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Figure 18 illustrates the effect of Lock number Y on the 
stability envelop for rotor with no preconing and no hinge offset. 
Unlike the forcing response, decreasing Lock number reduces the 
instability region. It can also be seen that the solution becomes 
unstable on the part of the instability envelop for Lock numbers 
of 6 and 3. 

Figure 19 presents the effect of structural daraping of the 

flutter boundary. With the inclusion of small structural damping 

coefficients of =^, = .005 for centrally hinged rotor v/ith no 

p 9 

preconing, the instability region is reduced drastically. Also 

for ^=^,=.01, this rotor becomes completely stable. This can be 
p 9 

seen from sim.ple linear solution (Eq. (106) that the structural 
damping coefficient has a very strong effect on the comparatively 
low direct damping term for lagging motion. 

In Fig. 20, the penetration lines of increasing rotational 
speed Q into instability envelop for different are plotted. 
Figure 21a represents the limit cycle flutter solution for cOo=0.8. 

P 

Points A and B here correspond to the boundary points A and B on 
instability envelop (Fig. 20) . One gets large fluttering oscilla- 
tions at the very initiation of flutter and then the amplitude 
reduces with increasing rotational speed. The solution becomes 
unstable for part of the speed range. It can also be seen that 
the flapping response amplitude is lower t^an the lagging ampli- 
tude mainly because of, comparatively high aerodynamic damping 
in flapping motion. In Fig. 21b, the limit cycle lagging ampli- 
tudes are plotted for different at the lov;er flutter speed 
side. The solution shows nonlinear stiffening effects for 

_ P 

of .6 and .8 and softening effects for higher Wq of 1.0, 1.2 

p 

and 1.4. 
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SECTION 9 

EFFECT OF VELOCITY GRADIENT 


In the earlier analysis, it is assumed that the incoming v;ind 
velocity is uniform and perpendicular to the plane of rotation. 

In reality, wind is not uniform but has velocity gradient because 
of earth's boundary layer. Generally, the velocity profile near 
the surface of the earth is approximated by a pov/er law relation C®, is] 

- = (^f 

l/o 'Ho/ (116) 

where 

V = wind velocity at height h 
= v/ind velocity at height 

p = constant guantity depends on topographical conditions of 
the place, between .167 to .40 approximately. 




Putting the above relation for the blade spanwise position 
r with azimuthal angle \p 


J/ 



- h. cos ^ 
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(117) 



is zero at the lov/er end position of the blade. 

Using Taylor's series and keeping terms up to second order, 
one can expand Eq. (117) as 




(118) 


where represents the velocity at the hub 


/h V 

' 4 /} 


The induced velocity at the blade is calculated by equating 
elemental thrusts obtained from momentum theory and blade element 
theory. For small built in twist 8^ , the inflow ratio parameter 
X is approximated as 




(119) 


where X^, corresponds to the inflow ratio at the hub and 




(a = rotor solidity = 


nc 

ttR 


) 


Thus, the inflow at a particular station on tlie blade varies 
as the blade rotates, causing periodic variation in airloads. 
This acts mainly as forcing function in the flapping degree of 
motion. 

Substituting this new value of A in Eq. (85) to obtain the 
inflow velocity at the rotor and modifying accordingly the 
resultant flow velocities U and v in Eqs. 87 and 88, and once 
again, on performing the various integrations for generalized 
aerodynamic forces one gets the new versions of Eqs. 92 and 93 
representing the equations of motion for flapping-lagging rotor 
with sheared flow. 
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Flapping Equation : 


+zp(^ +^c|5) + ^(i + e +->'p ) 


I* { Q ^ Q C-J + ^ Cg 4 C^ 4 |2 ^ ^0 4 ^ 

-f-<p^^c,^ +<p^^q^ -h$cj>^c,s-^<pcB, +(p\} 

X C05<k| ^ 2>7 + ^ ^8 ^ 

F Az -i- 4 ^<pl >,3 + $<p"' I >,5 } 


(120) 


Lagging Equation; 


00 oZ 


<^ <p - 2 ^(^ + 6 cp) + <p(e + 3 ^/) 4 0 ^ 4 2 <^ c^. 

^ (l~iS'-i<p^-*-S^c^) +GCo^H-(cp~^cl,^-i$p') 

= - X I'b^ -h/3^8^ 4^6^ 4^<^>g^ 4 4C^B^ 

4 ^ 4 + I p" g,^ B,^+F^q + ^BC , } 

cos 4- ^ 4 +^<P E^ £g ^ fs^ +^^cp 


— X rnc: (1/ 


S 


I ^^E„ 4 p Cf? £,^ ^ I E,3 4 ^ 4 c? cj>" J 


(121) 


The various constants in the above equations are defined in 
Appendix IV. 
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The important difference between these equations and those- 
representing uniform inflov/ case, Eqs. 92 and 93, is that here 
in Eqs. 120 and 121, periodic aerodynamic terms are present 
because of velocity gradient (the D and E coefficients terms 
above) . 

The nonlinear response of flapping-lagging rotor in the 
presence of gravity and aerodynamic forces with sheared flow is 
obtained by following the similar steps as applied for uniform 
inflow case. On substituting general limit cycle solution (Eq. 94) 
into nonlinear Eqs. 120 and 121, and balancing out various terms 
one gets modified versions of R (Eq. 95) , v/hich now consists of 
four separate parts. 

(a) contribution from inertia and stiffness, RI 

(b) Contribution from gravity forces, RG 

(c) Contribution from direct aerodynamic forces, RA 

(d) Contribution from periodic aerodynamic forces, RV 
such that 

R.J = Rr(i,}) + R^(‘, i) 

^ « ( 122 ) 

=■ RIO (<^) + RGO(t) J RflO(c) + X RVO(i) 

RI and RG are same as given in Appendix II for uniform inflow 
case and RA gets modified by replacing old C's and D's in Appendix II 
with new C's and D's of sheared flow. The RV which appears because 
of inflow variation at the blade is expressed in Appendix V. 

The stability analysis here is also quite similar to that of 
earlier case except that the periodic aerodynamic terms have to 
be treated like gravity terms. 
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9.1 Numerical Results 


Numerical calculations have been made for some specific 
configurations of Case IV to demonstrate the effect of sheared 
flow on the dynamic response of the rotor. 

Case IV. 


-.12 

o 

II 

w 

-e- 

Up = . 693 

^DW = -0035 

e = 

.075 

10.6 

X = .1 

a = .023 

= -.14 

®2 ^ 

.105 

.005 

Ca = -005 
<j> 

Cd„ = .012 
o 

a =6.0 

0 = 

0 


This case approximately simulates the characteristics of NASA 
Plumbrook 100 KW windmill blade with the assumption of complete 
rigidity in the torsional degree of freedom. 

The forced response results are plotted in Fig 22. 


Configuration A; 


Configuration B: 


Configuration C: 


Represents the condition of the rotor being 
excited only by gravity forces with uniform 
inflow. 

(.e. no velocity gradient, p=0) 

Represents the rotor being excited by only 
aerodyncimic forces because of velocity gradient, 
p=.167, R/h^=.625. 

(i.e. no gravity forces, 

Represents the condition^ of rotor under combined 
action of gravity and aerodynamic forces. 

(p=.167, R/h^=.625, (j.j^^^=.00 35) 


In Fig. 22(a), one can find that lagging amplitude is little 
effected by the velocity gradient and it is infact not possible to 
differentiate the results of configurations A and C. Fig. 22(b) 
shows that away from the resonance condition, the flapping amplitude 
is predominantly excited by the velocity gradient, and near the 
resonance condition the increased amplitude is mainly due to coupling 
•from the gravity forces. 

It should be noted that the presence of any tower shadov/ effect 
v/ould combine with the velocity gradient effect here to produce a 
greater excitation of the flapping amplitude. 
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PART C: INCLUSION OF FEATHERING DEGREE OF FREEDOM 

SECTION 10 

RESPONSE OF FLAPPIMG-LAGGIMG-FEATHERING ROTOR 

In the earlier analysis, it is assumed that the rotor has 
only two degrees of motion, i.e., flapping and lagging. In the 
present chapter, the third degree of motion, feathering, normal 
to the axis of rotation is also considered thus making the rotor 
a three degree of freedom system. 

The fomulation of general equations of motion is very 
similar to that obtained earlier for flapping-lagging rotor. 

The inertia part is already derived as a general case in part A 
(Eqs. 29-31) and the aerodynamic forces are obtained after modify- 
ing the flow com.ponents of part B to. include the feathering motion 0. 

The motion of the blade as angular velocities about the three 
axes C/ ? from Eq. 10 are, 

• • , 

COg =. XI (sin 0 stn cb + cos 0 sin p cos <|j) Q cos p cos 4> — ^ sinc|> 

' ^ (123) 

CO,^=i^(5inQ coscji- ces0 sin.c|>) — 0 cosp si'nc|> _ ^ cos4> 

♦ • 

OJ^ = Xi, cos 0 cosp — Q sin p <p 

Modifying Equations 87 and 88 for resultant flow velocities U 
and v 

U = h. (_a COS0 cosp + 0 sinp) R (- ^05 0 S(np sinc^D + sin© cosc^^) 

+ exi(siu 0 sinp siiacj) 4- Cos G coscjj) 

(l2 

l> = _ n h.( sm 0 coscj3 — cos <9 sin (S sinc|>) + ?ip coscj^ + cosp smcj) 

+ 2 SIR COS0 cosp - en SiW Q cos p 4 - ^ |n(sih0 s\nc^ t- cos 0 strip coscji) 

- p sinc|3 + 6 cosp cos cp j 

The A is inflow velocity parameter which for sheared flow is given 
in Eq. 119 as 

= \ + \ ^cosr (125) 


and n is the distance of 75% chord position from aerodynamic center. 



The generalized aerodynamic forces are 


^^+^£40 ~ — 

(^e)n«M. " sin^J'h.clM * sinpJ>idC 

(r is the same as ^ and integration is on the span 
of the blade) 


( 126 ) 


Again, taking the built-in incidence distribution {from 
Eq. 91) 

Qt = 0, + K±^ Q (127) 

h. + e 

Like earlier analysis, expanding trignometric terms to the 
third order and after performing various integrations and putting 
together inertia, aerodynamic and stiffness parts, v;e get the 
nonlinear equations of motion for a flappin-lagging-f Gathering 
rigid rotor with flexible springs at the hinge-point in the 
presence of sheared flow. 



Flapping Equation: 


?('-¥) + 04^ - + zfe - §y ~z 9 +04 

r^/fp 

= |: |f^ + Fi P + f| 0 + 1^ fr c|>% f| 0"+ ^Cf, + F^ po +. [| <| 3 © + Fc, p 4^0 -t- F,^ fe 

+ f^,t ’8 -‘- F,^B^* FiiP * fiM*?' + P|s® + fis P'f' * *17 PS + F;%^P + p;=l +8 +< 50 % 

+ Fj, 84 + ir^ gg + F^^pp*+ l^s p<() + 1 8 % 1^, 4>p4 +f^j| p’ 

+ - 1 - F 3 ^ 4 )e^-i. F 3 ,< 9^4 4 F 32 < 9 pe f F 33 Q 4 ® + ^ 2 .M 

+ F 37 1 Q 4^ + F 33 <P 4 F-^ G f^h tep + 

+ X C06V I l?o 4 I?, p t S,8 + Rj pV j5^ + R, p.8 + Rj 48 4 Rt) (148 

+ Rio p <9 -<■ f?li *4 (9 -f f?u 0^ + Ri 3 P + R|4 4> -^- R|S Q -<- PF* + ^'7 1^® '^^iS 

+ f?,q 4 >G + RjoG p + R„ 044 R^J 0 Q + R ^3 ^ ^ §, 

4 haT^p^ 4 R:^qc|ic|5* i + Ri\(9^4 +<^3a<9p® 

(12 s) 


Lagging Equation: 


oa oo a o oi o ^ \ ^ i ? 

4 >-0p -20^-S 4 -^^(p ~^Q(p.cp+Ry^Q)- 2 ^(p^-< 9 cp) -+cp(e 4 :i^ ) -)- 0 ^-M 9 c|> 
-^\/>42 (<^^c|)-2^cjb +GojSq'(cp-^cf>^~fS^^) 

= ~ ^ / Go G( p + G ^0 ^3 P ■*■ ■'■ ^s 0 '*’ p 4 ■*■ Gq p 0 G^ 4 -’Q 4 Gc^ G\q p 0 

+G„4^8 4G(2^Q ■''G(3'P'''Gm4‘‘'^Gi^Q ■*'G|^^4^ + Gi7 pQ 4Gi%4^p -+G,(^4’Q + G^o 0^ 

4 G^i 94 ) 4 G22® ® *®23 4 ■*■ ^;!S 

4 G^<^ 4* 4 4 G^qC^S + Gji 0 p4 4 G 32 0 pS 4 G 23 G pQ 4 G^i^ pc|) 4 
f G37 g> Q cp + 43 4 G^QQ 4 Gm p P 4 Gtj2 pS 4-Gi^^t|)Qp 

— X Co^) ^ So ■*■ S| P 4 $2 0 + S 3 -4- S 7 4* 4 Sg 8 ^ + P4> 4 S 7 p0 4 Sg <pQ -)- Sc^ j?>4^ Q 

^ S.oP*Q +S„cP ^6 4 S, 20 ^ 4 S, 3 ^ 4 S,^^ 4 S,S 0 - 4 S,(,P 4 > 4 S, 7 PS 4 S,s| 4 » 

4 45^0 0 p 4 S;j,Q 4) 4S,2 0G4S23PP4)4£,4 pp"4S^g\l(f4S^^p8'’ 

4 Si^cjp^) 4 S 2 gCpp^ 4 S:;.^cP 43 ^ 4 S 3 ^^e ''4 $3,9^434 Sj^Gpe 4 Sjj 0 c|)gJ 

(12.R) 
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Feathering Equation: 


-^cp) -)■ ^6 ^^ + ^2(p(licp - 4 - R„,( 9 ) 

— fi'J’ + *^l F( 4 Gj pQ F Fj4*® ■*■ (f^+ G»')f5't’S + 4 4?© 

(<^3 + G6 - i Fo) (f; - f Fo) 4?4 0^4 4 (F, 4 G,) P4.V (G^ - Go') 

4^5^© +fl3p4'^^t4'^P 4G(sSp4h504'F'^l5pP 4 F,l, 44 +^lfe'^ P4^ 

^6 ^4 •<• fi7 ^ 4''9 •*■ *^17 p p <9 4 f^g4p‘t’ + Gigi^P^-' f^c) 4 4>® 4 G,q<4 PG 

-*■ (^:iO +(^2,')^P<P ^1 (94"4G2j©J5© 4 4 ©3^14^ 

+ + ^is4®4 4^^3s4'©p4 iis4>4 4G2s40^+^33^ P *-G,^^ep 

4- X cos V' [So^ + ^o4 -<- s, ^^4 s^ 4 4s 4 ($3 - ^ So) p’ + C>?s- i R.O 4 S^) 

+ C^M /?6) <4*^ 4 c/?^ ~ ^(?o)4 ■'■ % ® j® + R^Q <p t- Sj +(Ri 4 Sg) p4® 

4 ^13 PP 4 s,^ 4^ + 44 4 s,g ©p 4 Sj^p p4 4 s,,, |3 4 r,^ P 4© 

"■ ^ (ISO) 


The various constants in the above equations are defined in 
Appendix VI. In the derivation of the above equations^ it is 
assumed that the aerodynamic center, the c.g. and the elastic 
axis lies at the same chordwise position and there is no variation 
of this position along the length of the blade. 
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10.1 Simple Solution 

Considering small perturbation solution about some static 
solution 

^ ~ ^ (131) 

+ = 1?, + $ 

Q ©0 0 


Substituting this in Egs. 128-130, retaining only the linear terms 
in the perturbations 6, 4> / 9/ one gets three linear coupled equa- 
tions. To simplify these equations further, one arbitrarily uncouples 
these equations and considers simple forced response of a rotor with 

zero twist distribution (0.=O) and with uniform inflow 

1 

Flapping Equation: 








+ - G sin cost)/ 


(132) 


Lagging Equation: 

$ ^ ^ ^ smtK-Gc^coS4' (132b) 


Feathering Equation: 


G * e { j V, f.. X I - II X I 


and static solutions are 


1+ e -f- p/ 

p 






So tj,* gp ^ X ( Qof, II ■) _ , - ^(1 4 |)<fe 


(132c) 


(133) 
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For ths forced oscillations assuining steady solution 


p ~ a, sin. -f b, cos 
^ ~ CLi sin -4- cos 4/ 

© - sia 4^ + bj cos If 


By matching and solving equations, one gets 


G(9,(e 4-y^) _ g/3 


li 



,Z 

(e + )^’/+ 1 






( e 3 f)r 


( 134 ) 



b 

3 ~ 




( 135 ) 
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From this simple analysis one can find that the flapping- 
response is effected appreciably by the feathering angle 0^. It 
is found that this simple solution gives reasonable estimate of 
the solution, but more than that, this solution is quite useful 
in giving an initial trial vector for nonlinear solution. 

10.2 Nonlinear Solution 

Assuming a limit cycle solution of Eqs. 128-1 JO for both 
forced oscillations as well as parametric resonance 

S is 3 a.sin.occt' -f h, casi>c<r 

^ ' (136) 

(p ^ -i- Sin oc If -t- h>2 cos PC ^ 

9 ^ (9^ + Q_3 S(Vi ocif -h h>j ^<^5 ocH' 

Here, a=l represents forced oscillations and a=^ gives 
parametric resonance. Substituting these expressions (136) into 
Eqs. 128-130 and balancing out the constant, sinonf^ and cosai/^ 
terms from each equation gives nine nonlinear algebraic equation; 



The various terms in the above matrices are expressed in Appendix VII. 

The solution of these nonlinear equations (137) is again obtained 
numerically by the Nev/ton-Raphson technique. 


84 



•10.3 Stability Check of Solution 

The stability check for the solution obtained through the 
Harmonic Balance Method is made as done for the flapping-lagging 
rotor by giving small perturbations to the steady solution and 
studying the growth rate of these perturbations with time (see 
details in Section 5.3). 

The pertubed solution is 

^ L^io + sin cc<f' -f. 

^ ^ 4 (^)] C^c + kM] cos oc^ 

9 = (9^0 +§c(^) + [^50 + + cosccr (13B) 

Since 6 , a, » b-« is the steady solution for which the 

CO 10 30 30 ^ t 

stability check is being made and 6^, a^, b^ are the time 

dependent small perturbations given to respective steady solution 
components . 

Putting Eq. (93) into basic governing Eqs. 128-130, keeping 
only linear terms in perturbations and their first order derivatives, 
filtering out the steady solution, once again on matching various 
terms one gets nine linear algebraic equations 



This is standard algebraic eigenvalue problem, the nature 
of the roots explains whether the solution is stable or not. If 
any one of the roots has positive real value, means perturbation 
will grow v;ith time and thus makes the solution unstable. 
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10.4 Numerical Results 


Numerical calculations for forced response are made for 
case IV (see Section 9.1) v;ith oig of 14.3 (i.e., Wg = 14.3 for 

two conditions of with and without sheared flow. In Fig. 23, 
the response amplitudes represening flapping, lagging and feather- 
ing degrees of freedom are plotted for different rotational speeds, 
for a rotor with uniform inflov;. Comparing with the results of 
same configuration with the torsion degree of freedom locked 
(i.e., u)g-»-<») from Figs. 22(a) and 22(b), one finds that there 
is relatively less effect on the lagging response but flapping 
response is appreciably effected by the inclusion of the feather- 
ing degree of freedom. The flapping amplitude is increased near 
as well as av/ay from the resonance condition and this is quite 
evident from the simple linear flapping equation (132a) where 
one finds that the forcing function for flapping is dependent on 
feathering angle. Also, even for this relatively torsionally 
stiff rotor (ojg=14.3) the feathering response is quite significant 
and it is of the same order of magnitude as the lagging response 
or flapping response. The resonance condition, however, takes 
place at almost the same rotational speed as for two degrees of 
freedom case. 

In Fig. 24, the response amplitudes are plotted for a rotor 
with sheared flov/ but with the same inflow A of 0.1 at the hub 
as the first case. Once again one finds that the velocity gradient 
has more prominent influence on the flapping response as compared 
to other two degrees of freedom. Away from the resonance condition, 
the flapping amplitude is generally increased and lagging and 
feathering response amplitudes remain nearly the same whereas near 
the resonance condition all the three response am.plitudes are 
reduced because of the sheared flow. 
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SECTION 11 


SELF-EXCITED FLUTTER RESPONSE OF FLAPPING- 
LxAGGING-FEATHERING ROTOR 

It is primarily intended to show here that the blade represent- 
ing the three degrees of motion, namely flapping, lagging and feather 
ing, in the absence of gravity forces and v/ith uniform inflow can get 
into self-sustained oscillations by the interaction of unsteady aero- 
dynamic forces with structural vibrations. The equations of motion 
for flutter of this system are the same as Eqs. 128-130 except that 
all periodic terms are absent (because G=0, X^=0, • 

11.1 Linear Analysis 

Linear flutter analysis is worked out by assuming that the 
blade response consists of small perturbation motion (B/ 0) 

about some possibly large static positions (B^/ 0^) 

^ = d ^ 

^ ^ t ^ (140) 

9 ^ Qo + G 


Substituting this in Eqs. 129-131 and keeping only linear terms and 
derivatives and after filtering out the static part, the character- 
istic equation of motion can be put in standard spring-mass-damper 
form 
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— 0 



9 0 




0 

Mas 


$ 

+ 

Q 2 ^ 2 ?. 





0 0 





rn„ 


Q 


*-^2 ^33 


0 

- 
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^ _ 


- — 









4- 

kzi ^Z2 kz 


<P 




^32 '’’'^33 


/*w 

0 




L 




(/m) 
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and the 

static 

equation 

is 





H„ 

H| 2 ^ 

'H.. 




— . — 

B + X-f 


H^i 

Hz 2 

H2 3 



— 

^ ^ ~ 1 " 


I 

j: 

*^32 

H 33 


0 c 




The coefficients of the above matrices are given as 

M„ = /.o 

^3 = ^ 

^2t= ^ 

1-0 

^^3. = t 
M32. ~ ^ 

M 33 ^ ^4 

~ ~ ^ ( ^5 ^ ^O ^ ^ *■ ^o ) 

^21 ~ -^{-(^13 -^ii^((, *■ 

^zz- ^icpHp ^ ^*^(8 •*■ 

^23 ~ ^ {^15 ^ ^ ■'■ 0Q ^22) 

q,= 
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— (+ e 4- 




k,^=. c^^-X(F,+^5oFs+^o 

A.. = 0o + |^CG, + ^^^3 + 

/?22= e 4- -+- ^^6 '^ <^oGs) 


^^3= ^ + XCG, + -J(9oQs 

^31 = + ‘^ 0 ^ 2 ) 

o 


^32= I - (^ + ^ ^ "*■ 

^33 = F.(/- 



^5 ^ ^ 0 *^ 5 } 


H„ = I + € -t- ^ (F'l ^ F'a ^ ^ ^ ‘^o 

Hiz ~ a ^0 ~ ^ *9o fs ) 

'^|5 = ■k'fe-i-C^^ *^-'"5') 

Hji = ifl. H-XfG, +IS 5 + ISIG,) 

Hji = e + + 1^ ('■6 &,, * § Gj ■uSt.G?') 

= i? + 1^ (G,+ 6!, S5). 


H 3 I = 

^33 - 

Hb 3 = 


i t - I- Go + S G| * t f; + ‘SVGj-' G-,) * g’CS3 - ^ )} 

fPo- 

<^>_|" + /?„(,+ Vg")-X{ ?G, + <6li + ?t(S;+Gs) ^■t%*‘k%’'s +f«oGs} 


Coefficients F^'s and G^'s 

\ 


are given in Appendix VI 
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Equation (141) can be rewritten as 


[m] 

[c]" 

c 

^ f 


O 

ro] 

;[M] 

-r 


s . - - / 


-- 0 


(143) 


where (M ) , [C] and [K] are respectively inertia, damping and 

stiffness matrices, each of order 3x3 and q is column matrix 


M = 


0 


(144) 


writing 


= KJ 


-1 

e 


(145) 


Substituting into Eq. (143) results in an algebraic eigenvalue 
problem which can be solved numerically by using any one of the 
standard subroutines. The solution will give six eigenvalues, in 
fact, three pair of complex conjugate roots. 


^ + {, Aj 


(146) 


gives the frequency of the rotor in the presence of airflow 
for the nth mode and represents the effective damping of the 

nth mode. The damping coefficient of the rotor for any mode in 
the presence of airflow can be obtained from Eq. (146) 


^ ( 147 ) 

The mode shape for each one of these natural frequencies can be 
had from the corresponding lower half of the eigenvector. Also, 
one can obtain the phase plot for these three degrees of motion 
from the Argand diagram of complex eigenvector. 

It should be noted that c=0 defines the critical flutter 
condition. C positive gives damped oscillations and ^ negative 
results in unstable oscillations. 
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11.2 Nonlinear Flutter Analysis 

Again assuming the limit cycle flutter solution of nonlinear 
equations 128-130 (v/ith G=0, X^=X2=0) of the form 


P ~ ^ -4- Q, sin b, cosocif* 

4* ~ sin DC 4^ 4 - bj COS c=c 

Q ^ + Q3 + b^ cos £xg» 


(148) 


where a represents the ratio of flutter frequency to rotational 
frequency. 

Once again, trying to find the limit cycle flutter solution 
for a knovm lagging amplitude by setting b 2=0 and a 2 as knovm 
quantity. Placing these expressions (143) into Eqs. (128-130) 
and balancing out the constant, sinai); and cosat/; terms and putting 
appropriately one can get 




— — 


— — 

'’u V ^ A6 




^o 

^21 ^22 ^23 . • 


a. 


A.O 

Aj, A 22 - ' - - 




A30 

1 

t 

1 

1 

t 

1 

1 


i 


\o 

- - - - - - . - 


z 


A So 

fl., - . . - . - 


PC 


Ago 

1 

1 

1 

1 

1 

1 

r 


0c 


Aro 

- . . . _ 



■ 

A So 

Si ’ . _ _ . _ 




Ago 


(149) 


Here Z = (w^/fi) 


2 


and the other terms in the above matrices are 
defined in Appendix VIII. 

The solution of these nonlinear algebraic equations (149) is 
obtained numerically by I'7ewton-Raphson technique - 
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11 .3 Stability Check of Solution 




The stability check of the nonlinear flutter solution is made 
by giving small perturbations to the steady solution 


^ +■ pcif +r^,o+ kMjcosocr 

A p- A . — n • (150) 

4 ^ ^ ^ b^(^)j cascc^ 

S = Sci,* 8t(v) + J/33c ■♦■ SjClf)] SKlCCy + [1>3„ + 

where b^^,, a^g, 0 ^^, 33 , 3 , and b^g = 0 with 

corresponding Z and a is the steady solution and B , a, , b, , (}> , 

Cj-xC 

a .21 fc> 2 / <i>^, a^/ b^ with same Z and a are the time dependent pertur- 
bations . 

Placing (150) into Eqs .• (128-130) , linearizing and subtracting 
the steady solution, again on balancing the constant, sina!|; and 
cosa<j; terms one gets 




— ~\ 


0 

Ti Ti2 Ts “ ■ ‘ ^ 1*3 


^ 1 ' 
ft 

T.. T22 - - • • 


S, 

t> 

T3, - - - • . 


t, 

0 

\ ■ - - - 


i 

1 

1 

' 


» 

t 

. 


0 






(151) 


Again the nature of the eigenvalues will decide whether the 
solution is stable or not. If any one of the root has positive 
real part shows that the solution is unstable. 
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11.4 Numerical Results 

In Fig. 25, the total damping coefficient ? is plotted for 

different rotational speeds for Case IV v;ith Wq= 6.0 (i.e., u„ = 6.0 w.). 

y u 9 

The dotted lines correspond to the zero structural damping configura- 
tions where as the full lines represent the configurations v;ith struc-' 
tural damping coefficients of ^ „=Q.=Ca~ • 005 . For any particular 

p <p y 

rotational speed, one gets three damping coefficients corresponding 
to three different vibration modes, but in the diagram only two 
branches are shown since the third mode is comparatively highly 
damped. In Fig. 26, the relative response amplitude for these two 
branches for the damped case are plotted. From these two figures, 
one can find that the high frequency branch I which primarily repre- 
sents feathering mode with small coupling from the lagging, is 
unstable from the very beginning for zero structural damping, and 
becomes completely stable with the inclusion of small amount of 
structural damping. Here, the lagging response is roughly 180° 
out of phase with the feathering amplitude. Low frequency branch 
II which is predominantly feathering-flapping branch with compara- 
tively small coupling from lagging is not much effected by the 
small amount of structural damping except that the instability 
boundary shifts to slightly higher rotational speed. This instability 
is like classical bending-torsion flutter, and it is also found that 
at the critical flutter condition, the flapping response lags the 
feathering amplitude by about 90°. 

Fig. 27 represents the damping coefficient ^ of Case IV with 

more torsionally stiff rotor i.e., w-=14.3. Here one finds that 

y 

the high frequency feathering branch I does not become completely 
stable with the addition of small structural damping but only the 
instability region shifts to higher rotational speed. Low frequency 
branch II which is now evenly coupled between flapping, lagging 
and feathering modes becomes completely stable with the addition of 
small amounts of structural damping. 

Fig. 28 shows the nonlinear limit cycle flutter amplitudes for 
flapping, lagging and feathering motions for Case IV with Wg=6.0. 

The bending over of the response amplitude curves towards decreasing 
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rotational speeds depict a typical nonlinear softening spring type 
system. For — less than about 0.3/ solutions becomes unstable. 

These nonlinear solutions indicate the possibility of sustained 
limit cycle flutter oscillations occurring v;ell belov/ the critical 
speed predicted by linear theory. These might be initiated by a 
sufficiently large finite distrubance/ and as such may be potentially 
dangerous. 



SECTION 12 


CONCLUSIONS AND SUGGESTED FURTHER WORK 

In Part A, the nonlinear response of tv;o-degree flapping- 
lagging rotor in the absence of aerodynamic forces is obtained. 
The analysis is made for forced oscillations as well as para- 
metric resonance under the action of periodic gravity forces. 

The effects of various parameters liKe initial coning angle, 
flapping to lagging frequencies ratio, hinge offset, etc. on the 
nonlinear response amplitude is discussed. 

In Part B, the analysis is extended to include the effect 
of aerodynamic forces on the nonlinear response of the flapping- 
lagging rotor. The inclusion of aerodynamic forces produce quite 
significant effects on some of the configurations for both forced 
oscillations as well as parametric resonance. Also, the self- 
excited flutter solution for flapping-lagging rotor is obtained 
after neglecting the gravity forces. The effect of the various 
parameters associated with aerodynamic forces like inflow ratio. 
Lock number, initial coning angle, etc. on the forced response 
as v;ell as critical flutter boundary is discussed. The effect 
of wind shear on the forced response of flapping-lagging rotor 
is also investigated and it is seen that the velocity gradient 
produces little effect on the lagging response but it has 
appreciable influence on the flapping amplitude. 

In Part C, the third degree of motion, feathering, is con- 
sidered thus making the rotor a three degree of freedom system. 
First, the forced response of flapping-lagging-feathering rotor 
under gravitational field and v;ith sheared flow is studied. It 
is seen that even for relatively torsionally stiff rotor, the 
flapping amplitude is very much increased with the inclusion of 
third degree of freedom and also the feathering response is 
appreciable. Then the self-excited flutter solution of this 
three-degree rotor is investigated and it is found that the 
feathering degree of motion is very important for the flutter 
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analysis. The effect of rotational speed on the linear as well 
as nonlinear flutter solutions is discussed. 

One may further extend these analyses to include 

1. Tower motion 

2. Blade flexibility 

3. Interaction of forced response and flutter response 

4. Subharmonic and Superharmonic response 

5. Improved Aerodynamic Model 

The inclusion of tower flexibility would involve the inter- 
action of the overall vibration levels at the tower hub with the 
supporting structure which may give rise to instabilities similar 
to the ground resonance problems in helicopter rotors. The inclu 
sion of blade flexibility would bring a more realistic model of 
long, thin windmill blades into the picture. The interaction of 
forced response and self-excited flutter response makes the 
analysis a more involved one and this interaction can become 
particularly significant v/hen flutter frequency is near tlie 
resonance frequency. Because of low structural damping, the 
subharmonic and superharmonic response may be important. The 
variable inflow model taking effects of tower shadow, and yawed 
flow, etc. and also the inclusion of returning wake and other 
unsteady aerodynamic effects will expand the domain of the 
dynamic problems of the wind turbine rotor-. 

The investigation of these problems would contribute to a 
better and more realistic understanding of the aeroelastic 
behavior of thin rotating windmill blade. 
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APPENDIX I 


The various constants used in the definition of aerodynamic 
forces are given below. 
n 

C, = Sin ©o 

R i 

Q = ^ sin 0o clh. 

o 

C 3 = sin 0 , 4h. 

(J - I + 

~ " - C2( As + - z 6. + CB2(A - ^S) + -^ CB^( A^S -I- ^ A - f?®) 

Q = -y - « 2 (Ae + 1 -)- i C 3 (aas <• S.) - CB, e icB,(A- h |e) ^-cb 3 ( As +A | - ie) 

g = -2 q(A - | 8 ) - 5 q ( A’e + A I - ^‘s) - CB, - cb2(2As + |) + cs^ a ( a - 5 |.o) 

C, = -C8,(I-^)-CB„AS - I 

C, = 20 ,( 1 -!“) t lAC^e ■*■ 2Cj|(l-|*) +CB,0 -C^A(l-|‘) 4 <0B^|S. 

C,„= C82A 
Q, = i CB, 

" *'B| ■“ T CBi |- 
C|3 — - CB, — ?A 

C„ = -c, +icB 2 A 

Qs= -Ci-’C^I <• icB^A 

Q, = c,o-e“)*jCjA(0-|.8^)4 sq |(i-s') 4CjA'6’'4 3CjA|.(s-|a^) 

4CB,(8-|a®)-CBjA(l-2e“) +3C8^ |.(s-|.g*)-cs aV 8 -§-S')+ 0 ^ A|(3s’'-l)4C|g(9-|g^ 
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R 

( r COS So 

1 

L COS ^ 

<1 

J*^h. <a, 
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R' 

M 

t^' o 


^ 

R? 


Cclo 

a 


(i+ 'i|')+^C.Q-Ci(X-|8')-2C3(X0 


, e.’-0'X -6 Cj,XQ 


+ 6CjX(A- 2-|S) 

, „ . , e e 4- ), s. - §■» ^ s'* ' ^ 

- -M ^^{\^ + f A - 3 ■» V j 

+ 5 BC,6 - ^ BCi(A- ^©') 

r \(0 4 -C e.(i_^') 4- ^ BC.6 - 

^ = C,(j- ^ ) V XG -k- Ci 3. } 

a"- % \s + ^ ^ ^ -6*)- ■'■ ^ B®") 

-Qo(R-e + %XQ'^-3s SR ^ 

’<1 "a ^ ^ 

^ =i - Cj,X - ^ 

\0 

= .iC, ^ BC,A 

B,^ = -^C-c, I + 8C,X 

r? _ ^ \ A- C| 

D,3 3 a 

R - - - -2 CxA 

m CL '‘. 

\,;...^.c4«.»'v.t<.-»-'>i-'.i‘'‘-t*'>-‘»''--'';f-‘‘-»’'i 


CB, = 

B, + 

Cdo 



a 



iL C^lo 
S a. 

CB 3 = 

63 4 

a. Qo 

a 

6C, = 

B, - 

-L CJo 
^ 0 . 

BC:, = 

B;.- 

. 3-- CJo 
3 -5T 

BCj = 

B3- 

- CJo 
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APPENDIX II 


The elements of response matrix 'R' consist of three separate 
parts, i.e. , 

(a) contribution from inertia and stiffness, RI 

(b) contribution from gravity forces, RG 

(c) contribution from aerodynamic forces, RA 

so that 

= RI(i,j) + RG(i,j) + I RA(i,j) . 

= RlO(i) .+ RGO(i) + ^ RAO(i) 

The various elements are given as 

RI(I,I) = 2 + I + e - I - 

= 0 

KrC'.3')=o ^ 

= e? - Nj { + 

RIO, 5) = 0 
Rl0,6)= 0 
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•Rl(2,0=O 




•RI(2,i)= WpZ +O-«:"')0- 
R1(2,S') = ~2cc cJpJz (p 


RI(^iM) =0 

RI(2,S)= 0 - 2tcl - Nj { 

RI (2,4')= - 20 c (I + 


Rrfs,o = o 


2cC C0^-f2 




RI(3,^') = o 

RI(3,5')^2ccf^ + Q<<2) 

T?I(5,0= 8- - Nl{«('^>^.-'^*'“''>} 

'RI(^,^') =1 0 - Nj I a,q^ + b, bj.| 






Rl(M,S) = o 


Rr(‘<,4')=o 
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i?r( 5 ,i) = 0 


SI(s,2')= = + 

■RI(S,3)= 2cc(^ + 0c})') 

■RKs^')^© 

RI(5;5)=: Z-ve +Q"-a“" + Nr 

Rl(S^6) = -2cx:y2 

= 0 

RI(6,1') = -2oc(|+ Qc^') 

W,i)= 9 - 5 |<tj_ 

Rr(c,M') -0 

Rl(i>,5)= :ioCjJ ^ 

RI(i,i) = H+£ -K= + e^-g^ +Nj|A«=(3qf+l,';).i(af+jt>^^ 

•RI(!?(|) = -c3p Z Nj.{ia,Q^(i+oc*)j 
T^IO(l') =: O 

■RIO (3) =. 0 

Rlo(^):= -Z c|3 - Nj. I o:(9 (a,b^- a^b,)| 

Rro(s) = o 


RI 0 C 6 ) = 0 


Gravity terms for forced oscillations 


»?G( ),1 ) = o 

Rg(i, 3'^= GCi- ^)- Ni I 

RG(i;S-)= 

RG(i,4)= Nj[Ga,^,| 

RG ( 2 , 1 ") = — I ” 

RG( 2 ,^):z 0 
RG(2j3) - 0 
RG(h^)^ - Nr { 

RG(i,s) = 0 
■RG(?;4)= 0 

t^g(5;|^= g (/- i*’- _ iVj. I ^(a?--v 3b^ ta^ 

RG C3,^) = 0 
RG(3, 3)=:0 

RG(3,M)= -Nj I ^(a.Q^tSb.V^ol 

RG( 3 ,S)=.o 


RG(3/) = O 
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RGCM,6)=: G(-L_p^ .i 

_____ _ 

RGCs.i) = , w r ^ , 

^ b,)| 

( S/-5 ) =: 0 
RG (s^ 3 ) S5 0 
Rs(s,4)= 

^GCs.sr)^ 0 

'?&C^;4)=0 


+ 








<f6)-^) = 0 


R<^(G,3)=i 0 


RG^6,i^) = G(/-|^_ 


R<3a,S)cr cJ 

R^7 ) = Q 


/"g + 3b^)j 


1 no 



-RGOn) = ° 


1 ^ 60 ( 2 .) 
RGO (s') = 
T?60(M) = 
1^GO(S) : 

RG O ) 


-G0f(- i") + ^i{ -^3a^ + b^)} 

' 6 z z l. 8 

0 

G(f- f- - '^l{ ^ ' ^^(3Q,« 

: - Wi I ^ Qi (a, bj, + <^Z^’,')| 



Gravity terms for Parametric Resor.ance 


1^0(1, 0= N, ||(aNb" 

■RG(l, 2)= O 
•RG(l,3)= 0 

T?G0,5-)=o 
t') - 0 

RG (•?, 2") - £) * 3“?)] 

•RG(2,3)=^ 

•RG(2,i/) = 0 
RG(2,s')= jGf 
RG(z,(.)zz {Q<j’6 

= 0 

RG(3,r)= iG|0 

RG<l,i)= &0:- £- £) - r,,( * 3i>‘ ■)} 

-RG(3,^i)= 0 

Rg(3,s)= 


RS(3,0= -iGgi 
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'RG(M,i)= 0 
= 0 

T^G(^;b)= mJ |(a"-fcj" + cN b^)| 

RG(b,S)=o 

RG(4/6)=o 

T?G(5,l') = 0 
liG(S,2)^ T^f ^ 

T?G(S,3')= 

■RG(s,‘/')= o 

RG(S, 5>= -sff - l’- I') + ^ ^ )| 

1?G(5,4)= -XG(cJ- sg) 

RG(C,l')^ 0 
■R&(6,2) = iGcjbe 

RG ( 6 , ^ ^ = 0 

^G(<-,s)= -ke(i- a^) 

Rsu,o= nj-f- i) <■(>:')[ 
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>?Go(0 = Nj.{ ^(a,b, 4- 


■RGO(i)= O 
■RGO(s) = o 

-/^r{ 

"RGOCS)::^ 0 

■Rgo ((>)= 0 


In the above elements 

f 0 for linear case 


N 


I 


1 for nonlinear case 

. 
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Aerodynamic contribution terms 


= -N^ |ec(a,b^- 


■RR(t,i)= o 


T5fl(l,S') = o 


R/^(i,i/)= - Cc^b,')(C»3.- ^i)} 


Rfl((,S) = 0 


■Rfl0>6)= o 


RAC?,!) = o 


R<\(1,3)= ccC, +ec|t^C„ * f f C3al+b*)C^ 

*ffoUb:)c„} 

Rflf2,H) = 0 

RHfts)= -VN„{ f a,b,(c„- c,,)} 

Rft(i,4)= «.c, * £c^(j>^c„ + ocg'c,,, +N*{-^(<iNb'')C,„ + ^a,ajC ,3 + £(;a‘,.b’jf 

+ £Cft^b:-)c,3 


(?A(3,1) -O 

RA (3,2) = - ore, - »C|<5 C„ - «|3f„ . _ w^|s(a,‘-t;)C„ -*f (<i;t3i;)C,i,.£(a.\-t;)C,jj 

R()(3,3)= _,^C, 

Rfl(-;,i() = o 

RA(5,3)=-»cc, - - <C|V„, - S(«?-b;)C„ A- f l,liC,3 +sfa,'*rf)C„ 

HH(M)= -%[ (C,,- C„,)| 
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'Rfi(^,^) = 0 
■Ra(V,5) = o 


{ cc(a,b^~ <Z,b^)(8,^- la')! 

^A(i|, s)=o 


RA (b,6) =o 


Rft(S,\)=zO 

Rn(S,i)=-oce^. N„{-f(oM:!')B„- 

(s,Lj) = 0 

Rn(s,s)= 2.5 s^ * |s^ ^’^n{-i«‘^.k(B.,- S,,")} 

/?/) cs,0 = - »" e, - or<iie,s *'^«{ f ('o'-bf )e„- f 'S.'^z B,, - £(31!f + fcf ) g,^ 


F.f\(G,t) = o 
Rfl (4,W)~ 0 



/?fl('6,S)- oc 6^ +cc^^a,3 -f oc^* 5 ,^ -t-ccf^ + 8,3 ^ S(q?- + 3b^ ) B 

+ + 1\')B,3} 

Rft(ti6)= 5<|Z z-g E, 4- a,b, Qa,^- B,^')| 
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1!ftO(,')= -Co-p^C;- t'q- -'^4 

-i^c"cQ,a^4-b.b^')CB, +icc"(al + b"^c,^ 

T^A0(2) = 0 
RAO(3') = o 

+ -i «c (<1(^3, + b, C, + -b cX. (Q-i. + b^^ — -b £>C + b^ ^ BC,^ 

RAO(s) = o 
RAO(4^=o 


In the above elements 



0 for linear aerodynamic case 

1 for nonlinear aerodynamic case 



APPENDIX III 

The elements A. . of Eq.- (113) are divided into two parts 

X 3 

(a) contribution from inertia and stiffness, AI 

(b) contribution from aerodynamic forces, AA 

Thus 

A.. = AI(i,j) + J AA(i,j) 

xj o 

A^q = AlO(i) + I AAO(i) 

The various elements are defined below 


AI(»,3^=<5 

Ar(/,5}= 

Ar(/,6)= 

BlCz^i) = 0 

fn(z, 2 ) = 1+ e - e"- n^{ +i.(af+b^ 

ftr(z,3) =. -zee • 

s) = cj^a, 
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filfs.l) = 0 


)- \ +e~oc~ +-^(a^-v 

f)I(3,ii)=0 

/?rC3,s)= W^b, 

/)I(3,4)= 0 

0 - Nj{a,a^| 

= 0 

fll(%3)= 0 

e +8^- g^-N,|i(i-oc*)(af+b;')J 
fll('s,l')= 0 

flr(S,i)=r 8_a|^ 

» 

fll(s,5) = 2 ccf^ + 0 ^) 

RKSjii) = 0 
pi (5,5)= Qj, 


flr(5,«)= 0 
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= 0 


9I((,,3) = 0- 

fn<c,^) = 0 
ni(c,s) = 0 
m(c, ^)= 2.J^ 

fno(i) = o 

f)io(z)^ea^- 

fuo(3) = accftj, 4 <9 <j^) 

flro(^) ~ o 

/?rofs)= Q,f - 

fnO((>)~ ~^if£<l,b,<Z2_{i + oc^)j 


In the above elements 


N 



f O for linear case 
1 for nonlinear case 
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Aerodynamic contribution terms 


fifttt, \) = o 

RR(i,2) = 

+ ^ih') + ^ t ^ 12 . %^)} 

M (l,‘^) = 0 
=0 

RO(I,(.)-0 

0 = -QaC:^ 

RR(Z,Z) = -2gCs-^d-7 

flnfj , 3 )= K(Cj •* c ,;)\ . 

RR(Z,5)~ 0 

RR(z,C)- I '2 C^io ■" 

flfl(3,l) ~0 

nfl ( 5 , J ) = -< c ( C , + +£<^ 1-^01 

fl« 3 ,J)= -JJCj -.|C, 
ftA(3,4) = 0 

RR(^S)=o 

ftf\ ( i , C ^ - - 0.2_(Cf^ ■*■ f ^13 ^ 
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+ ocQ,')] 

flfl(^,3) = fo bJ + I'O} 

f!RCM,s) = 0 
Bfl(M, 4 ) = 0 

f)fl (S,[') =z Qa 67 
AA(S,2') = 2 | + <^6^ 

ms,z) = B,„ 4. |‘S„ + 4 O - n 4 f fa'+il )6„ * 

flfifs,!-,) = aa^e/ 

Rn(5,5) = 0 

flflfM) =cc(6, 4.|‘e„ 4-^X)4-N*jS«;(^^4^j') 4.«(« 

flflf{,3) = 6j 4 . 6, 

M( 6,‘l) = o 
m((.,s) = o 

fllHi.i)- a, (6, 4 g B„ 4 |‘g,^ 4 ) 4 Nfl { i 6,s 4 ^ t; ) 4 
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f ®/3)} 




Q3.(<lT43b;)8 1 
4 ’ 



f^no(i) = - (c, ^cplc^ ^ ±a\c^ ^ Cyclic,] 


ftRO(z) ~ o 


flfld(3) = 0 


flflOC.)= B. *|‘b^ + 8,-. _«\af^b^) +(afH-b;)|s 


RRO(s) - o 


flf\0(C) = o 


In the above elements 

0 for linear aerodynamic case 
\Ji for nonlinear aerodynamic case 


N, 


= < 


121 


APPENDIX IV 


The various constants used in the definition of aerodynamic 
forces with shear flow are given below: . 

Q, = Co- A,(CB, * ICB.-) 5A,9(CH, + IcO-EB^e -► e, eS »A,e\c6, + IcsOt 

-JA,8 *(ch, + |c,) 

Cs = C5+ A,(cs^<- i|C8,')-S(A,CH,--^EB^) 
g = Q -E A, ( I C6, + i CB,) -A,6(S CH, + 3 |C,) + EB3 S 
C,= C, -sA,(CM,t |c,)<-E63-JE,9-sA,8(CB^+|CB,) 
g = C3-A, cB^s 

t, = C, A,(-CB^ * 1 CH,S V ^ CB^e") 

f;. = C,o + A, ce^ 
c = c 
c — c 

^13 ~ ^(3 ~ ^ 

Q4 ~ ^14 "a ■^1 

^15 “ ^(5 


B„-^A,(CM, + |C,)-EC3 tU, |oQ(l + |)-^E,0 +5A,0C6C,4-|BC,) 
^ec3)- .a,Q^Ccm,. 1 1) - I 

-JA, 8 ^(bc,^|bc;) 
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B3= Bs-A,(ch, H-EC3-i£,e-A,S{8Cc*^|BC,+ |^-). 

5= jtA,c,-i^A,8(| + |)-e(iE,+A,BcO 

g,= + |)-(e, + =>A,BC,)-Ja«f,8 +aA,8(ch,+ 

Bj= Eg +A,(-2Bc^-hecH, + e*Bc^) 

1,= B, •+ \,(CH, * I a.6 
B,„= B,.-A,ch, 

B„ = B„ + A, BC^ 

8,. = B,, + A, BC, 

B,.= S„-| %A, 

= B,^ - i A, CH, 


Ij,= -AjCcb, + J «i') + ? Ai8(C, + EB^8 +Ej6’ + |CBJ 

-iA,9^(c,+ |cO- 

Us = A,(cb, + l I cB^') ->. 9(-AgC, + i EB,') 

II = U>(CB, +l|CBj-A^0(2C, + 3|q) 1-EB^e 
I), = -2A^(c, + |c,)+EB,-2E^8-AA;,e(cB,*|cs,) 

= -A^'^B,S 

= Aj(-CE, + 5 C,8 + i; t8, s’) 
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70l(I> 



^,0 = 

P„ = 0 

D = o 

IZ 

J^3 = 

= ^b, 

l>s = i CB, 


Eo ''JCi * |c4-Ec, + + k I) + E^g ^ jA,S(Bc,+ I gc^) 

+ e’C % E, * EC,) - 5 A,s\c, * I CO- i a. A, e^(, * k |) _ | ^^ 


-|A,e"(BC, + |BCO 

Es - -A^fc, + 1 1 Cj) + EC,- ^ E,0 -A^8(i -t ec, + 2 I gcO 
- H ^.<^2 - % A,G ( H- II ) -8 f i E, A, sc,) 

II) - (e,-2^ A,SC,)-3^ ..?A,S(C,+ 

^ = A,(- 2 bc, + qe + 8c, s^) 


~ (c, + 2 . 9o Q ~ 1 c, Q* ) 

^to ~ 


i sc, 

= 4 sc, 

E = - 2 9lo 

d 2. 
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t-14 “ 


- i C, 


^IS ~ ® 


CH, = 


CH, = 


BH, = 


BHi 

CB^ 


K o 

R 

^6 J sifi B, dh. 

^sj h.'" COS.I9, clh. 

So 


■R 

J h= cos B„ 

K o 


BH, + y. ^ 
' 5 a 


BCu = 8H, - I- £ds 


~ R* (r) j =''>8o =l*|- 

? I ^ ^ ®° '='’^ 


fB, = E, + a 

t = E 4 ^o E, 

T 2 CL 6 

CCj = £, - i £s 

EC, = - £ C^ 
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APPENDIX V 


Element Matrix RV for FORCED RESPONSE 
RVO.i) = o 

RV(I,Z') = -H% P<.'C„*45>,0 - '^0 + 

RV(i,3') = 

■RV(i,H) = o 
RV 0> s) 

RV(i , s) 

RvfJ.O = + ■*•^(“1 

RV(2,2) = 0 

RV(2, 3) ~ 0 

RV(2,^) = - \{^('Q.a,-b,tJ(l)3 + 2l>J + + 

RV {2, s') =0 ^ 

RV (^>4) =0 

RV(3,0 = + i z(3^zb, 

RV(^,^) = 0 
RV(3,3) = 0 

RV(i,H) = - N^{^C^.b,3),^+ i- 
RV(3,S) = 0 

RV(5,6) = o 


= Pc t{(^,a, + 5b.bOD,3 

*K^-b>)D„-b;r>„-Ab,b,i„] 

= - ID 7 


- 14 ) 11 , 5 } 

- a, b^') j 
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RV(^,I) = o 


i ( Es I4i E,o ^ ^ 4 E,0 + i n 4 (a.a, 4- 2 b.bj E„ -V K+ b, )E, 

+ (Q 2 .+ 3bj') B,2 - b E,j - 5 b, b^ 

R\/(^.3) = ?c ^s -k t bq 
RV(%M) = ° 

RV (M,5) = i(E, + St, e,3 <- |‘ E,^ + 4 E,s) N^-E{ 

*(«?+3b^)E,^ + (tt: + k’)E,5-4E|c.-=>=,'^LM 

1?V(4,6) = t ^6 * -5 ^ S 

RV(S,I) = Nj{J;(a,Qj-b,b,)(E^+5 E„) •*-i('lJ-bL)E,3} 

RV(S,2) = 0 
RV(S,3') =0 

RV (5,M) =: ^a{ T r(^'^i-'^k')(Ei3 + ^ ^'0 -*^t) 

RV (5,5) =0 ^ 

RV (5,4) = 0 


Rv(t,0 = N«{^Ci«.b,-0^k,)E„ +tt,b, £„ +iO,^E,^+ 

RV(&ii) =0 
RV(()3) = 0 

RV(i,h) = Ns{Ea,b, E,„ 4 -L(5<»,b^-a^b,') E,^ 4X(3Q,b,-a,b,')E,5 + 04b,E,5} 
RV(4, S ) = 0 


RV(4,,6) = o 
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T^VO(l) = 0 


R[/0{z) 

Rl/0f3) =. -(n>^iA 

^\/0(^) =. 0 

Tiyois) = N^[ 

KV0(6') = £o -(- ^^£5 -f- ^A i { ) ^s '*'(^^ ^G 

+ («, + 3b, b^)E^ j 
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Elements of Matrix RV for 


P^ARAT»iETRlC RES0^7A^Jrn 


RI/(i>2) 

Rl^O/3) = 

RV 0,5) 

RV (1,6) 


= O 


+ ■!!>„) 


= 0 




=z 0 


RV(2,I) = 
RV(2, 2) = 
RV(2j 3) = 
RV(2,C,) ^ 
RV (h S) = 

RV 0,e) = 


+iin, 

2- Ms X 


Rv(B,2) 
RV(3,2,) 
RV(s,^) 

RV (3, s) := 

RV (3,c) = 


= 0 




0 
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= f^A{ + 2 E^') + 3; Si} 

RV(^,2) = O 
RV(^,3) == o 

= N^ { ^ Si) 4 ^'o 

R V (M, 5) = 0 
RV (^) 6 ) = o 


RV/(S , I') = 0 

RV ( 5 ; 2) = -(^ ^5 i i ^t) 
RV (S,3) = i(6s * S'® 


RV (s,q) = 0 
RV (5,5)= 




•RV( 6 ,i) = 


RV( 6 , 2) = 




RV ( 6 , 3) = 
RV ( 6 ,M) = 

RV (e, 5) = 
R V f 6 , 6 ) = 


= g £5 t 


(£<1 *f£,4 + #£,s) + '^Ai;{‘’>A£.3'£'’' ^'■<■"'’>■^'5} 
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RU0(2) = 0 
Rl/0(3) = o 
RVO (^) = - Nfy 


•RVO (5) = 0 
RV0((>) = 0 
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APPENDIX VI 


The various constants used in the definition of aerodynamic 
forces for flapping- lagging-feathering rotor are given below 








h = 


^ = 


f 




^5 = 


hi = - 


2 fA.c, +A,ch.) +2 f(A„Cj+ A,c,) *06, -(a;c6j+ EBj) + 2 |cb^+ ce, 

-<^| -|q +(■'1/8.* -I, CB,) ♦ i t(A.CB3 + A,CB,) 

-|<^2-f.C3 + i|(A„CB3 2-A,CB,) 


F - - F 


c, + (a!c 3 -V E,) - 2 |-q- ^ C3 +2(AoCb^ + A,cb^) + i |-(XoC63v a, cb,) 
2 (A„q + A,CH,) -2|(A„C3-f A,c,)-cb, +(A^cb3+ EB3)- |cb^ 


F 7 - 


F. = 


F,= 


Fio = 


fT. = 




-I^AoCBjH- A. CB,") 

- 1 (C 8 ..|C 63 ) 

-2(A'cj -►£,■) *2|.q-n ^‘Cj-2^(AI.cb^4A,ce,)-2 |(A„cE3 + A,cs,') 

- (A,C2 A, CH,) - i C8, 4 ■!( A‘c8j+ EBj) - I CB,^- i SjcSj 

-(A„q +A, cH,)-2 |(A.Cj+A,C,)-ics, 4 iCA'cs^ 4 ES3)-|cS3-J.|>3-i|2 

-i (A.C, 4 A, CM,) - i |(A.C3+ A,C,) - f CB, + 1 { a'cBj + ES^) - J f ce^ - J 


^3 - ~ (CB, + & C^) 

^ |^ 2 -(tCB ^4 A,CB„) 
Fis = -I (c% 4 I C 83 ) 


F,t= |(|CB3 4CE3) 
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F„=-lCB. 

fiq = 2 ( A 0 C 2 . + A, CH,) -*- CB, + |- CB^ 

cr 2 e c, -h ( Ai^B^^-SC, 

*20 “ R 

fiz = -i(AoC 83 + A,C6,') 

= (Ao^^2+ a, 

F,, = ^ CB, 

Pis = +f cB^-i p.s 

Fa, = 

= -2(A,q-^A,CH.)-'^®. 


Fas = -^' ^ 
F2, = 


« 7 i 

0 

II 

— r — S- c + 

F3, = 

~ 2 ( Ao Cj^ + A, 

II 

(M 

_C8.-|C8, 

% = 

- (AoC8, + A, 

^SS = 

-C8, 

F35 - 

-± CB, 
R ^ 

Fss - 

CB, 

^37 "■ 

-CB, 

Fss = 


It 

<r 

UI" 

-ifis 

p.. = 

-2C, 


\ (Ao 

CH,') 

C8,) 


F = --^ F 
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G = 9 o( i + i §- + 2 W ( Ao Q + A, CH.) -1- K + A,c,') -( A„ BC^ + eQ 

^0 a. \ 3 R J ^ 

G, = i { I Cj - 2 ^ "• *^4 

G. = % ( I A. * I A, ) + ^ I ' t 

+ 2 |-(AoSC 3 4 - A, Be,') 

^ —2 

63 = - &(l* i|)-('AoC.*A,C<,,)-i|(A,C3^-J,C,)+('C8C3+EC3)-iBC3 

Gs = - ^f/+- J |- + 2 :|^') 4--^ (^Ao -*- Es')"^(Ao^2+ A,Ch,) -2 |-(AoC3 + A,C.) - 

+ ( Ao 6C3 4- EC3) -Z^BC^-^^ 6C3 

= - ^ (l-A. + I A,") -i" f -(A’.C5 *E,)*|Q-^(A,bc,+ A,bc,') 

G,= i(eCj+A,C,) *’|(6C,4.|BC3) 

G, = i(Cj*|C3)-2i(A.6C3+A,BC,) 

G, = -2^C2A^E5)+f 


4- 2 BC, 4- 2 f 6c^ 


s,„ = -i ^ ( I Ao * I A,) it - -k ( Ao Cj t.) + |Cz * i fl G5 - CAo BC, + A, sc,) 

-2 f (A„bc3 + A,Bc,) 

G„ = - f ^( 1 A„+ I A,) -2 ^ |-(2A„+A,) + iC|- 5 {A„C 3 + E,) *g-Eo.* i 
-(A„6Q+A, BC,) 

6„ = -f ^(|Ao - I X,> % I % M " %C,-|( A‘oC 3 + E,) * IIC3 

+ |S? C3 -i.(A„6c, + X,BC,)_ i. a (A^BCj + A,sc.) 

= C, 4- I - ^ ( Ao^Cz + A, BC^) 

G,4= 2^(1+ + ( A0C2 + A, CH,) 

G,s = |(C 2 + |C3)-=^|rAoBC3.A,sc,) 

1 -^A 



G,^ = (XoCz + A. CH.) -V 2 BC, + 2 I - i C^ 

^>S “ ^ ^ 2 . 

Gr = ^ ( t ■*■ I - I 
^21 = ^ I <^z “ ^ ^ 

*^22 = ^(AoC3+A|C,')+2i. ( SCj. + ^ 

R 

G,3 = _(A,q4-A,CH.)-2BC, 

G,^ =-iG ^ (A,6q + A,BC^) 

G^5 = --^C, - ^ Cj; + ( AoBq+ a, BC^') 

I'^O ( AoBCj, + A, BC^) 

^27 = - ^ (|Ao + |A.) + C, 

^26 ^ 

r = - H- CJo €. 

^2C) 3 sr R 

G30 = -| % I -•k(^o=X--^'‘^-* 0 -'%° 

S3, = |C. 

S33 = -^(l''"-*- 

Gj 5 (4C^ + A, CH,') + 2 6C, + 2 I sq 

= c, 

S33 =1^3. 

Gj, = - C, 

^37 = 

r 

CL 
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G.„ = o 


= - 1 6C. 




Gm - - 2 1 BC, 
= - 2 i 6C, 
= -2 8C, 

^HS — ° 

~ ~ 2^/ 


^0 - KC^^i •*■ ^ ^Bj') 

/?, = 0 

2i,c, + 2|A^q-EB^ 

Rj= A^cb, + x I A^cb^ 

= ^ A, I CB, . 

Rs = 4-5A^C8, \^cb^ 

R 

i?t = -?A^c,-a|X^C^4.EB^ 

R, = -| A^ CB^ 

Rs = ° 

= -2(£, 4- A,C6, 4-| >^CB^) 

R|0 “ ~ "2 

Ru -2A^C, -S|A^C^ + E8,^ 

R,a = -^A,c,-i|A,c, + i EB, 
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Ril^ ~ 

/?I5 = 0 
1 ?,^ = o 

Rfj = - CB, 

Ri^ = o 
Ri<^ = 

Rao = A^CB, 

Rj, = o 

R 25 . = - Aa^Bj ^ 

R 

R*2 3 ~ A2^B, 

R 2 M = ° 

Ras =° 

Rzt = o 

Ra*f ~ — ijA^Cj 

i A^ C6, 

Ra<i ~ a Aj. C1B| 

R 30 ~ "a Aj C6| 

R3, = --?A^c, 

= o 

R 33 = -Ai C8, 

So = Aj C, + ^Aj^Cj^— EC;^ 

s, =-2iAjBq 

R 

= 2 ^ A^( 14- ^ I) ^ + U, BC, + 2 I A^ 
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£, = - i,C, - i I + 




EC. 


Si- ^ Ej-5 A^c, -2 £. AjCj^ + EC., 

5c,= - I % |ix- E^-2 A^®c, 


^7 " ^ -^1^2 

Sj = -2 I A 2 6q 

Ss = +24 C,h-2| X 2 C 2 

S,o “ - % A2 - 2 ^ X^BC ,- 4 £2 - X2BC, 

Sn = - ^X2-| % f A2-CE2-X28C, 

% = -i«.A2-i<^A|A2'|E2-%A2(eq SC2^) 

S |3 = - 2 A 26 C 1 

5,^ = A 2 C 1 

S ,5 = -7 I A, 6 C 2 


S,J, = 2 IA 26 C 2 
^(7 “ ^2 C, 


S ,6 = ° 


5„ = 

S^o = — ^2 ^1 
5z, = -2A^8C, 


'ZTL 


= 1 A,c, 


Sj _ — Aj.^, 


— ^2. 8^1 


Sj5 =r A2 Be, 
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BC, 

^2S " 


S30 “ — "k ^2^1 


S3, - o 


c = — 2 A-, 

^32 a 


^23 ~ ^z ^1 
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APPENDIX VII 


The elements of response matrix *R* for flapping-lagging- 
feathering rotor consist of four separate parts 

(a) contribution from inertia and stiffness, RI 

(b) contribution from gravity forces, RG 

(c) contribution from aerodynamic forces, RA 

(d) contribution from periodic aerodynamic forces, RV 

such that 

R^j = RI(i,j) + G RG(ij) + I RA(i,j) + J RV(i,j) 

R^q = RlO(i) + G RGO(i) + ^ RAO(i) + ^ RVO(i) 

The various elements are given as 

Rr(i,z) = o 
Rl0,2)=- o 

f 

<9c - |(a, oc(Q:, ~ ^3 1 ^')! 

0 

RI(l,q=^ o 
Rl0,7)=^ 

RX (i.-i)- 0 
Rl(l,<i)=^0 
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Rr(2,0 


KI(Z,2) = oj^i -N,{i(Qi+b>;) + 

+ L ec\al-hl) 4 shl) +i')] 

RI(Z,3) = - Zocl3^.{z - Nj: {c<a, 0.^1 
Rl(Z,H)= ° 


Rr( 2 , 5 )= 9-2^^- Nj|b,b^(,+oc’)} 

RI C2,6) =-2PC('yS+e<^}-N,(cC£Z^a3} 

Rr('2,7)=0 

Ar('2,2)= <^('/- <9c - ^i{h^>k0~^ 

RT(2,q)~ Nj.^{(K(3a'^-hb^) -i. ± PC^aJ _ ')| 


/?r(3, 1) = 0 

^I(:3,2) = ZCCOJ^^ 


/?r('3,3) = oj^i + (i-c^"-)(i-4)-^t~^- + ^ - N^{i(<^f+^<)} 

4- -*- 3b:)~i ) - i- 4 4 3b;)} 


RI(3,H) = 0 

Rl(Z,S) - Zoc (^ -h Sc) 

RI(2,Q)=^ 0 - /Vjf-^a,a2,('/ + <x^)} 


Rr(3,7)=o 

RI(3>,S)= -2ocf-N^{^(CL^t-3l>^)- ^(a^-b^)} 
RT(3,^) = i(l-od^) ~2 ^Gc -Ni^^(l-oc^)a,Q^'^ 
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= 0 ^ - Nj I (a, Qi 4 b, b^) pcCa.b^-Qib^)! 

RI(Li,Z) = 0 


(?r(^(, 3 ) = 0 


e + Oc - - N£(- 5 (l-Pc')C^r+ ) - 'z(^s -*- *^3)(l )} 


R1 (tt,S) = o 
Rr (h,e) = 0 
RI (^n) = Alr|(Qxa3 + 




RI(^,^) = ° 

* 


RT(5,{) = 0 


Rl(5,Z) = {2 ^»^zC'-+ 

/?rC 5 , 3 :) = 2 ccf^ 

Rl( 5 ,Li)= 0 


RT(5,5)= Z + e -«: 




4 ^( 3 Ql + b. 


.=)} 


/?r('5,6) = -zocjz - /Yr{«ca,c?3j 


/?I(S,7) = 0 

/? r f S, 8 ) = 4 pc^ ) + 6)^ + A/i (i 0 + ^ 3 } 

Rl(s,8) = 2oc(Rm9 + ■*- A/f { jcc(’^^i^ 2 + 
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RI(6, I) =0 

Rz(G,2) = Nj(^«:(a^a3+ 2b^b3')J 

RI((>>3)= 2|(^ _/v^{ia,a^(l + cx’‘)| 

RI (i,^) = 0 

Rl (^, s') = 2 oc +• fC b, 

Rl (R,^) - ^ +e -cc 4- 0^-^ -i- CC (30^ 4 b^)-X 4. zlf) 

-TfOC (3a2 4-b^) +^(«3 +3b>3)| 

RT((,,7) = o 

Ri(^, 2) ^ N,^±cc(a,a^-^^kK')} 

RI (^,9) = ^Ji-hcc^) 4 -z4iG, 


■RI(7, l)— - iVj I Cj^ + b, b^l 
RI(7, 2) = 0 
RI(7,3) =0 

RI(7,5) = 0 

RI(7,(>) = o 

RI(7, 7 ) = >?„ (/ + 2 O)^ ) - A/j + 6^* ) - bi^)j 
RT(7,S) = 0 
RT{7,9):^ 0 
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= 0 


RI(Z,2) = ^ 

RI(%,3)= -2oc-b^^n,{l^ (a\ -yl^_ lcc(a^^ 
Rr ( 9 , q ) = O 
Rl(2,s^=. i(l+cc'^ 

Rl (s,(,) = -2 oc(R^q^ ■^^<5^) ~ N^(<^a, Oil 
Rl(s,i) = 0 


RU^>^) =R^0 cc " -to ^ ^ 40- 

Ri(8,q)=. -Zcc^^fs R^r;^ 


^")~N,{^(l-i-oc^)(st + 

(^l~4)} 



Rl(<^,l) - o 

RI (% 2) = 2oct I loc K- y ) - i cc ('a,V b : )} 

Rlf % 3 ) = i-Nj_f iQ , Q 3 ( c < C)J 

Rl(‘i,R)=.o 

RI(%5):= 2oc(R^i9,-i-^<p^)^N,{pcb,b,] 
RI(%7 )=o 




+■1^- C'-«’)(a^-J- ) _ ^ ) -1 cc^qI _ 
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)} 



fjro(0 = 

RlOfa) = 0 




0.1 ■+ 




^10(3^) = 0 

rio(M = 


RI0(5) = o 


RlOr^) = 0 


rI0(7)= 


cc")(a,a^-^- k h^) + 


RI0(S)= 0 

RJO(^) =0 


>^cc(a2bj'-Q3b,')j 
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Gravity terms for forced oscillations 


K&(i,i) = 0 

RG(i,^) = 

RG(fj5) = -f- 

RG(I,C) ^ 

RG(l,l) = o 

RQ(I,%) = +/f ^+ sal 4.^1 + 

RG(h^) =0 


RG( 2 ,I) 
RG (2,2) 


+^7.K) 

o 


RG(2,3) = 0 

RG(2,H) = 


RG (2,5) = o 
RQ(2,G) = O 


RG(2,l) - -( -^ 




i>: 


/?e (2,S) ~ O 


RG( 2 ,^) = ° 


«G( 3 , I) = 


RG (3,2) = O 


RG (3,3) = 0 

1^0 (h‘i) = 

RG (3, 5) = 0 

O = ° 

RQ(2,1) = 
RG('3,S) = 0 
RG (i,^) ^ 0 


RG(^,\) = o 

/?G(^,2) = 

R^eh^) = -iH 

RG(^,H) = 0 


RGC^I/S) 

/?G ('i/,6) 




RG (^, 7) ~ 0 

KG(4,?) = -i(S,-Rt) ■* 'V,{^(3a,.Jx + i.,fc,)| 

/JG('/,‘i) = N^{l(a,h^ + a,k,)j 
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RG(S,2) = o 
RG(5,3} =- 0 

RG(s,u) = N^(^(3a,a^ + b,b^)-^(Q,b, +a,h^)j 
RG(S, S) = o 
RG (SjC) = (^ 

= if +/v^{ i7C3a,a2-hb,b^)J 
RG (S, S) = o 
RQ (5;R) = 0 


RGC^.l) = ^{ -^(^zb3+a3b^V 7!-(a,Q,+ 3b,b,)} 

RG a,z) = O 


RG(G,Z) ~ 0 
RG (6,5) = o 



RG ( 6 / ^>) — o 


RG ( 6 , 7 )^ + 


RG (6,%) = 


RG (Q,^) - 0 
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RG{1, \ ) = 0 


RG(7,2) - -I ^ f (9c) 



4- 


i- 

16 




RG (7,3) = 0 


RG (7,^) = ° 

/JG(7,5) +N,(ifc,b^| 

RG(7,6) = 0 


RG(7,'7) = o 

RG(7,S) = 
RG (7,^) =0 


RG (S, I ) = 




RG C^/ 2) =0 
RG (g, 3 ) = o 

RG (g,^) = 


RG (g,5) = o 


RG(S, O =° 

f^G(S,7)= -f- 

RG (^/‘S) = o 
RG («,"t) = o 
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RG(^,2) = 0 

R&(^,3) = 0 

RG (<\, S)^ 0 
RG (q,6) = o 

RG(^,7)= N^{ ^A)] 

RG (S,9) =- 0 
RG 0,9) = 0 

RQO(i) -0 
RGO(z) = o 
RGO(Z) = 0 
RGO(>h) = o 

RGocs) = / - - Ni/ rC^Oj ) -» + ^3)} 

«G0C0 = 

I 

RGO O) = o 
RGO(S) = 

RGOO) - - Qs^O} 
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Gravity terms for Parametric Resonance 


R&{ 1 , l) — Nj I ^(P 3 b,.+ Qifc' 3 ') ~^i') ) I 

RG(t,z) - 0 
RG(IjB)= o 

RQ(hh)= N,{-!^-(Q,b^-ta,y,) ^■L(Q,Q^-b,b,)j 
RG(i,5)= 0 


R G (l/'j')— j if ( k + + ^3^3} J 

RQ(hi)- 0 
RG (h‘})= 0 


RG(Z,\)= 0 

R&(z,2) = ~i nd^+ii) +Nj. 

RG{Z,i {)=0 

RG(^> 5 -)^ 

RC(Z,(,)= 

RG (2i 7) = <5 
RG (2,S)~ 0 




RG(5,f) = 0 

/?g(- 3,3)= T-4-rr+#)-'v,{7ktNfC} 

RG (^,^)=^0 

RG( 2 ,S)=lcpj^^Nj{^h,b^] 

RG(i,G)= 


/?Gf 3 , 7 )= O 


RG(i,2)= --^ 
RG(3,R)= O. 


4 (V, [!■(•«,% i’) ■* 
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RG(‘i,2) = o 
RG(H/Z') = 0 

RG{HiH)= Nj[±(^l,fc34Q3fc,') + -zi^l'i)} 

RQ((<,S)=^o 
RG(H, G)-o 

R&Cii,?) = N^{^(a,b^+^2^,)j 

RG(‘i,S)=o 
RQ(^,^)= 0 

RG(S,l)=uo 

RQ(5,2)= 

/?GfS,3)= '^x{iM+k^s)} 

RG(S,H-)=0 

Rc (5,s)= + + 

/?G(S,6)= -i(>^-^<Sc) +N^{i^a,(i3l 
RG(5,7)~ 0 
RQ (5,2)=° 

RG (5,<f) = -t(^c-^<Pc)->- ^ ' 

RG((nl) = 0 

RG((,,Z) = 

RG(6,3)= ~^ip^ 

RG (^i‘i)= 0 

/? G r 6 , 5 ) = - X ^ _ I ^ X 6, f>3 J 

RG ((>,&)= { €) “ 72 ^ } 

RQ (bil) ~ 0 

RG (Q, 9)= - ^(<S'c~ ^iktj 

RG(i,^)= 0 
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RG(7,Z) = 0 
RG(7,3)=o 
RQ(7,H)~ Nxllfaifc, 

RG(7,5)~o 

R& (7,6) ~0 

RG(7,7)^N,{i^(a,b^+Q^b,)] 

RG(7,S) = o 

RG f7;9) = 0 


RG(S,l) = 0 
RG (i,2) =0 

/? c, rs, 3) = - i * i f g’ + + s? ) + ^ i 

R& (ijR)=.0 
RG (8,5) = o 
RG (9,6) = “ 

/?G (Sf7) = o 
RG («iS)=o 


RG(^.l) = 0 

RG(<i,Z) = - ^ + '^^c) ) 




RG(<),3) =0 
RG (%‘i) = o 

RG ('>,5)= C<9c-|<^) + 'Vr^^b,b,| 

RG (^.6)= 0 
RG (^,7) = ° 

/?G ('‘f.S) = if(90 + 

RG (9,7) = o 
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RG0(1) = o 

RGO('2) =0 


RG,0(2') - ° 

RQO(H) = ~^r{ 7 f ^3^3')} 

RG0('5) = 0 
RQO((>) =■ O 

RQ0(^)- -Nj.|i(Qzb3 + ^3^2')| 

RQo rs) = o 
RG0(‘i) = ^ 
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Aerodynamic contribution terms 


RA(1,I) F,^)~ 


+ -i- 




RA(i,2) =z 0 
RA C, 3 )=z 0 


'^Af'.‘<) =-^{ i(a,a, *hh^)F^ 

+ f (»3 fc, - «, fcs) Aj, + f ’ (a.5. + t, t-n 'i7 


+ +C*^iQ 2.-^ k,b^') f^s j 


RAO, 5 ) = o 
f^A (lj 6 )-o 


= - (^2 + t f^c-^ f^, ■*■ <^c F,^) - N^{i(^,a^■^■b,b^)F^^ F(aJ^b])Fo 

+ J K + )hl ^1 ( 0-2 + bj )F,^+ PC (a, ) ^4 + «■' r 

- f(a 2 b,-a,b^)F,^ ^f(a,b^-a^b,)F,^] 

RA(i,S) = 0 
RA ( 1 ,^) = o 


RA(2,l) = o 

^ +S.F, t *Sy|£K7o + <^'/5s) +F^fo 

* T« t.bj/Jj +^ccY 0 i 03+ * ^(35j + tl) F^sJ 

M f2,3) = «(F,3 j. <|b Fs Fr t <53+ I"/;, -t «A {f 

* f *C)5s *f f -f('«i-fc])53j 

Rfl C^/^) =o 
RR(2,5)^-(3<P,F; 

* i '57 - 4- t bj A„ * 4 *>.‘>= } 

M(Z,e) = "^('Fs <■« F»*«:F, *;?l§F7*®"'i8*€'i3 + ®‘ F30) " ^ 'is ‘‘f ‘'•“‘'jS 
- ¥ «. <I> F3, - E (3=,S t;' ■) F„- S (-<31 3 ) F, - bj -jllo * f(»3'- ‘1 153 * f ‘’=‘>3 fl} 

RR(2,n) ~ o 

RR(2,S) = -('^->-^^F^ + ^F^^(pfy^^(pf^ + ^^F^ 4- + 3(^" /^j) - b,b^/^ ^ i>^)Fa 

fif) (2,R) = ocffTs ^22 Q.,0.^ 

«2«3 Fio - b,b,)F^, -ioca.a^F^:, -f Q^a^ rj3 _ ^ + ^>z ) ^<? / 
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RA(-^, 0 = 0 

RA(3,2') = - Pc 4 f^as ^ ^ab) - '^a{ f ■*■ +<3. ~) 

^ ^(:al-.3b]')F^^-. b.b,F,,+£(«l-b])F3,\ ■ 

RA(b,3) =-('F, -^2p^r3^<^^4-(P,F;^64/^ ^2;6(?/7„-^c)^V^s)_A/4iff«aQ3^^^a03)^ 

+ f Q,Q3 ^ ^ ^ C«a + 3b* ) F^^j 

RA( 3 ,^)= o 

RA f 3 , s) = — oc ( Ft^ ^ F,^ 4 . Ffc^ 4 - F^^ + ^ 4 _<^ Ftj ->-i 3 c Psa) — ~^i ) ^a-j ~ ^ 

f b, f^a, S 3b0^8 + f bl ) f roj ^ 3b;)p3„^ f f«|-b])f^33 - ?b>3 P,,} 

RA(3.0= + ^ ^ ^3 fT, 

RA(-i,7) = 0 

RA( 3 ,s) =- cC^f;^ + < 9 c ^2 

- f b,b3F3, 4^(a,a,^3b,b,)F3,-^f b,b3F3a-*-^b,b3F33+^(Q;^3b:)F3,} 

RA(3,^)= -(Fa + 2^cFs + |^ + <ib/^+^cbb- 4 )*To 

4- A f 3b: ) ^ (Q3% bj ) /r^ _ f fO. ba 4 Qa b.) ^, - f V«' ~ ) ^^34 ~ f >:-b:).J,-sVz^.} 




=. G, + ■>- ^fi,ji('^2^3 4 bjb 3 ')G<, + (’Q/Q 3 + b,b3')G|o4 ~(a,i\- ^ab, )Gj3 + ^-C^jb, ~ Pf bO&aj 

+ Xe>c (Qjfcj- Qj{,^^G 3, +i»cY‘*i«3 4 b.bj'jGj^ 4 oC*(aj«3 4- b^ b3') G^^^ | 


RR(‘i, 2 ) = O 
RflC'</ 3 ')= o 


+ 


Nfl I lf‘<,'^3 4 4 (^^2^3 4 (\b3) G„ -4 «T('a, b^ - 0^ b,.) G^g 4- ^(Qab.-Q, 

f r«3^,-«.b3)^3, + f'KS 4 b,b3')G37 4 -f b,b3)G„3} 


5 )= O 


(^, 6 )= O 

R ft (un) = G, 4 G,„ 4 4 4- G,^ 4 -y, 

4|(«3*+ b:)G,^f ^(«.t^-«3b,)G 


{ i (^, ^2 4 b. ba ) Gq 4 «/' 4 b^ ) G,o 4 i-(Q/ 4 b' ) G„ 

2fe 4 <=^(^ 2 ^ 3 - Qjba) G30 4^(036,- a, ^)G22+ rKba-Qa% 3 j 


Rft(i^,S)= a 
RR (if,‘>)= 0 


\ 
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Rfl(5,0 = o 

R/}(S,2) = G, +J^G, + +■ ■*■ ^Si^cGi^ + /^a|-^(S 52Q3 "* -^- 5-b,b3G|0 

+ I OC%,b.^ G3t ^ ^ CC^(Q^Q^ ^ 3 b^b3)G4q - Gqs } 

Rfl(5^ 3) — Gj^ + (9^ G ,7 G^^ + Gjg +Gc + ^ Gq,"^ — W/ a,<52 ^23 

■••■^ r^/^+ ^/ K ("SQj + 4) ^3S ■'■ ^3 ) ^36 “ ^ (‘^2-bi ) ^iV~ Y ^1^2 *^3 8 

f - ^b)'^32 } 

RR(S,^) = o 

RA (S,5) = ^ ^ ^(5cS + 5 ^(9c G^I -f- %{ V i 

*- ^ (Q 1 G 2 + 3b,i^')Gq3 ~ 

RR(5,(>) = - Oc(Gh^ ^- ^cG,q -^/f ^ ^;27 ^ZR + <^3o) + ''^A ) ^23 

-^Kfl,a,G,s- f Q,Q^G„- ^(3Qf^6^)G^g - ^(aUbl)G,ci-^(3a^^k^)Gso 

^r(‘'i-‘‘i)ea,} 

f?Af5,7) = o 

/?fl('s,s) = G 2 + 2 ®,Gs+/fG 7 +^Gj 4 ^<$Gq +'v;,[^fc,b,Gq 4 ^(sa^bf )G,o 

+ 55-r3<?i + 4' )G/, ^ bj ) G,i + f r<?/4 «>.) ^3^ ^ f V^N ) ^36 + Q2 - ) ^37 , 

*■ ^ b,bi (Gqj + Gqq) 

Rfl(5jR) = -*-Pc^Zo + 4i ^31 + +/i<5^^3f ^31 ‘S 33 ) +^A{^Q.fl3<2r,6 

f C3a,Q2 + 6 ,b,)G 3 , - fa.Qa G 32 - f 02 ^ 3 ^ 33 } 


^53 

= 30 


T5fl(6,t) =0 , , -N r / , ^ 

RA(4,Z) = oc (G,3 + <6 G„ ^ ®cG,7 ^3S + 1%*, ^ ^ " i ^ 

-t ^C4"-^b’)G2q +^CQ2-*-3 4)Gis + f r«3 +3b*)G2t-.-eC«|-4)G,7-^b,b^G^g^S^a|-(^')G3^J 

RflC^.S) = Gj + 2^ Gj + ®c^7 + ®c Gq + ? /I 6J; G/o +A/fl{-q (Qi<^3 ^ 4 

+ ^^a,Q 3 G 3 , + + 4t^')Gqq-f^‘^zS<^M3} 

Rfl(^,R) =0 

RR(Q,S') ^ cc(G,^^ G,g + G,^ ^ G,7 ^ G„ + (j6*G„ + ^ ^a { ^ («N ) & 

4-S('«l-b* )G33} 

(?flC4,6) = ac^Gq + /f G^ +0cGg -f;f<9^Gq+a<^^®^G„ -/-/VAj-^Q.QsGq +^c<^Q3Gm 

+ CjQj ^3^ + fsa, <23 + 4 b3")Gq3 - ^^Q ,«3 Gqq ^ 

RA (4,7) = O 

RA r6,S) = ^rS.s + /|C^2o^<^G;2,+®cG,j4/|<^,G3,^,f<9,G3,+<?&5',G33')4..V,{-3:4fc3G,g-^b^4G 

■*■ ^fQ|Qz -^^^4)^3, 4 - ^ b.bj G 32 + ^ 4^ *^ 33 ! 

‘/?fl (6,9) = G^V2£iGs4.^G^+qbGs +^"g,o ^ ^®*G,2 + f ^^G,, +2f(«N34')G,o 

. 4 («,V34* )S„ . |-(«3* f ij ) G,z - Gj, - 4 )G3^- )G37 - ^ 
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e/) f 7, 1 )= _ Go - 6U + Gg) - + G^) ^^(6s- f«) - (<l^Gs - n 4 1 (a.a^ + b, ( P, 4 Gg ) 

+ ('fliQj + b,b3)&7 + (QpQz + b, +Gt~.^) + 3:(^l-*- b^')(Fi + ■*■ | C'^^’ + b* )(Gi - &) 

+ ^(Q| + b*')Gs + f (a.t.,-Qzb,)G,^ (a,b 3 -Q 3 b,')Gn +‘>^Kb,-fitb^')G,g +^(Qjb3-Q3^')G,q 

+ «: {'Q 3 b, - Qj bj ')^' 4 G, ,') + «c^Q 3 b, - Q, bj) G^o 4 ^^(a,Q^ 4 b, bj') ^ (a^a^ + b_^ b^) 

+ ^ (al 4 bj ) G3g + ^Y°3 ■*■ ^3 ~) ^<<0 4 4 b,b3'^ '*' b^ )^‘/6| 

/?fl('7,2) = O 

/?b(7/3) = 0 

/?fl('7;‘/)=— /b 4 ^3 4 (^7 ■•■^g} +(^2*^3 ''' 2 ^ 

4 ^ 4 b,’ ) r^-l 4 G^ - ^ ^ 4- I ('C' ,.blXfi,'~f)^i (^1 4 b] ) Fs 4 (Q,a^4-4bJ(F, 4G<,) 

4 x(a,b^-Qib,)f:^ +^(a,b^~ + K(’fl,b,-a,b^)f^g + ^C- 22 ^ -Q3b,)f^^ 

4 ^r«3b, -^|b3)Cf,o4 G^,) 4 «:(« 3 bj-Qzb 3 )/r^ ^ fV «>'^2 + ^b,)5,,4«\Qza3+ ^abj') ^35 

RflOiS) = o 

RR(->;6) = ° 

RR(7,7) G^ _<^Vg -A';, / 4 -(-Gg) 4 X(^«i^b;)G7 4 I r «2 4 b*) Fg +r 0 z« 3 ^b^t 3 ')|r 

+ {a,Q2+b,b^)Gs 4^(Q,b^-3zb,)fj-^ +^(Q,b,-a,b,')G,<,4|fa3b,-a,l^)G2j.4s(a3b^-c,t^')F^^| 
/?/l(7g) =0 
RR( 7 ,^)s =0 


t?A(«g,t') =0 

Rfl(9, i)=-&o-^^ Gp, ~<^F, 4 (P^Gj ~^/§^c^7 ~^H^0^3'*' '^6~-§-)~^(^(, 

- ^)~®cG5 — 4 b>ib3')(F^ 4 Gg') 4 :jb,bjG^4i b,b^(F^4 Gg_ ^ ^ 

4i(34^4b^)(F,4G^)4i(a4^b;)(G3-Go)4-if3a|4b])G5 4^H,b,G3,, 

+ «Vz-bt) 4 -^ 4^V<4 -^fK«3 4b^b,)G3^ 4«\4 43b^"^G3g45Y«jV3^")G„,-f" b.b3G,,} 

Rfl ( 8 , 3 ) = «: C ‘fe F 3 4 I G ,3 4 ^ <56 G,g 4 <p^fT^ ^■<p^(i, f-j + ^S^ G ,7 4 - A//, {-fTa, a^G/^ 

- s (a," + b*)G^, -f(3al-^bl)f^^-foc c,Q3 ^ «,Q3 g,, 4 ^ Fg 

+ f Q.CzG,g 4 ^ f«Nb3)G2l4 ^ Q.agGjo} 

RA(8,^) = 0 

RA ( 8 , 5 >-Fo -{| F, -S^F^-P.Oc (F 7 4 Gg:) _|Yp 3 + ‘^ 6 -%)- 3 t^<- 

^ f * 1 | (’^.4^14) ~ *'/a| 4 b,b3(Fi 4 Gg") 4 ^ b,b^F^ 4 -{^(scif + 1 ^ X^i + G^- 5 ) 

(F,'- S ) 4 ^ (sal 4 b; ) fT 4 i b,b, C»i 4 64)4 ) G3H f ^b. b, f^4 
" ^2 )^^-f'b,b3G33- ^*b.b,&3g} 

/?ACg,6) = ccf|G,4 4CjbF./ 4|<jbFs4^,"G,g4<^(9,F9 4^<aG^)-/V^/^C«f-b^)G,^4^Q,a,fi6 

4 K Q, Q 3 ) <5|6 - f ^22^3 f ^3 4 S <5, 
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^flf5s,7) = 0 


^ + ri +^^<?cC’s}- Wa{ 4; b,b, ( s + <r.s) 

^ -^(3^1^+ bf )Qq "m C 3^2^ 4 b^ ) 4 a b, 4 ^ Ca,b^+ *^:?\') 

4 ^\al -bl)C^^-^o.n,h^Q^s-f b. b3 G^, -. ^\af- bt ) G^2] 

4 ^ (ai^-bf ) G,7 4 ^ fa,^-b^ (q. q^_ b, b,) G, ( ^q.Q, 4 b, b,^, ( F,, 4 G^,~) 

- ^ a, a3 G, 5 - ^ (sa^ 4 b,^ ) Gio - ^ ( " f 4 


RA(q,\-)=^o 


RA (^.4) = - oc (<^ |T^ ^ (I G,3 4 4 4 4<^ (P, 4^ 4 Gn ) " { I b, b, G,6 

4 ^ b; ) G^, 4 S CG* 4 2b/) f;, 4 ^ c< b, b^ f7, 4 ^ b. G„ 4 ^ (a ^ _ b* ) Fis - ^ b, b^G.s 

- f ■’■ ^ C‘^3 ~ ^S? - f ^1 ^3 


4 2/|G, 4<^F; 4(9^G, 4t«Pc(F7+‘Sg') 4«/|<9 cCS7 4 2^<^(F3+G^- ^)^<jl((i^G^) 


+ 3/|V53-|s)4(P/g^| -/VA{i(‘^*‘^543b^b3)('F74G?) 4 1 Q, Oj G^ + ia,t32(Fj4Gg - & ) 


^ ( K 4 il>/)(fl^G^ ) 4 ^ (G,^ bf)(-G3 - Go) 4 XfG/4 3 b] )Gs 4 

4 fV‘^/4b;)Gv,4^V3G,G34b^b3)G3s4^V3S^ 


■^oc’a.Ga’^sif 

4 b*)G3S 4^(3-lf4b )G (^0 


Rfl(‘i,^)=0 / 

= -off/fS4c^fr^ 4^<^fT^4.,f"^,j4<5b(p, /79-^-f'? G,^)_ A/a/ 

" ^ ^1 ^ ^7 4 |F b,b^ F,^ 4 ^(Gi^4Sb*')Gi5 4- ^ bjb^ F,q 4 bib^Gjcj - ^ 

R Arq. 6 ) = - f Fo 4 g F, 4 5> F^ 4 g G> (F 7 + Gg) 4 ? t®c F 9 + gY Fj+ - 5 ) 4 ScfeYf"/- &) 4 «/ 4 ^§^(F,fGt^) 

3. ^ 

-WA{ija.Q3(S+^<&') + ■^lifC^i'+sbtxPi+Gt-S V ^(^iXb^)(F;- E) 

-}^(^|4 3b3)F_4 -La,Q,(F^ 4Gq)-^\aXb^ )G3,t^+ -^tfC.^a,G 2 ^ 3 ^+ '\s 

4 " 4 ^ ("Qa* 4 ) f'SS ~ "4'‘’^^^I*^3^3S ~ 2 ^ *^|Gi ^ 2 g | 

Rfi(q,7) = O 

/?AC9,S) = - cxr{ H^IS 4<$f;s 4)^^ ('/io+ ^ 2 /) +/?^C ^ 2 ^ 4-^^G^O 4<;^(Pc j 

-^n{^C^,ci^-b,b^)fr^-^f(ci,^~bJ')G,7 ^f(^/-b/)hc) 4 S ('G,^, - b,bj G,7 

4 ^ CQ,q, 4 3b, b, ) r^o 4 ^2f ) ^ ^ f r^/*4 3b,^ ) G^o 4^C«2'4 3 b|) b^, 4| b, b^ 
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= -{ + -^Gg) -f •+ ■+ ac^fij -+^^<5. Gg j - ^a{ ^ ^i^z(^r*-^i) 

■* ■— (Q^4- sbf ^ Gj -[f (Qj + 2. b^ ') Fg ^ {(X\h^-^ QjbO(^o + ^?l) 

- - 1 ; ) F3^ _ s V, a , G 35 - C3 G40 - } 

RflO fi) = - f Fo F 3 . <5bV^ 4 (P/^ ^ ^ ^^ 6 ) F, -4 S, S } - I ^ ( a,^ + bf ) F 3 

4 ^ (a’ 4 ti ) F^ 4- ^ (q| 4- b] ^ Fg 4 ^ (Q,q^ 4 b.b^^') Fg 4 -L (q^q^ + b, b^ )F^ +-L(Q,aj4 bb 3 )p 3 
4 f(a,b,-Q,bOF,t4f(a,b3-a^b.^fi-7 4f(a,b,-a,b,')frg4..^(Q.^-^3b0^.q 
4 I (a3b,_ a, |(a3b,-Q,b3^ F,, + ^ 

RRO(Z) = O 

Rf^0(2)=o 

Rf) 0(^)= ('Go + 4 (P/Gs 4 ^ G^ 4 y| ,p^ G 7 4 ^ Gg ) +/Va [ ^ bf ) &s 4 

+ jC^/ 4 bj* )<3g 4 -^(diQ^ 4 b,i^)Gg +- 3 (^i<^i + b,bj')Grj 4As)^8 ^ 

4 ^(Ril>^-'^3b,)G,'j + ^(^zb,-Q,b^)Gis 4 2('q, bj- <33 b,)G;ci 4 a.cq 3 b, -Q| b 3 ')G;?o 
+ «('«34-q,b3)e^, 4^Y«,Q;j-t-b,4)G3z, +q^‘Jz<^3->-kbs')^ss 4fH’-+b*)G3s 

4^V«I 4 4 f V^,33 +4i>3)G^2 ^ «Y< 4 4’)G^^ J 

RRO (S) = o 
RfiO(6) = 0 

f^^OM = ~(M"g, +P,<^F; ^pJ,G^^c^S^F^)-NA{±(a^4.b^)G, -^i((i,Q^^b,b,)F; 4.j(Q,Q^,b,i^)G^ 

+ i.(a,Q^^b^b^)F^ +f (a,b^- q,b,) f)3 4 1 (q^ b,_ G, b,) g,,, + f (a^b,- < j., b^) q,s 

4-|(a3b,-a,b3^F;3j 

=0 

RflO(9)~ 0 
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Periodic aerodynamic forces terms for forced oscillations 


RV/(i,0 = o 

RV(h2) = - 3 -*■ ®c ^.7 + ^ ^23 + ^ (Q.Q^- 2 b, bj ^23 

= 0 

/?\/ (l>5) = ~i (R,^ + ^Rli-bScRjq + ^(^R2'i+^ ^3 o)-^A S ^23 ^^^5 

+ ('G,«2 + 2 b, bj)^27 +3bj*')R29 + ("qI + b’ 4 (a| + 3 b’ - b,P3^3(-^^33^ 

T?V/ (/, O = - J ( ^ ^ ^6 + 4 ^ 4 4 l9c ^,|) - ^A 4 ( a,q3Rq + (iCi^as + b^) R„ } 

Rv (b7) =0 

RV (/,g)= ^§/?20 + <^Rai + 4'?*2+ itRsi +^®C^32+ <^^b^3 3)-^X 

- 5b,b3R3o +('^.Qa+3b,bjR3j + (Q,a3 4 2b,b3)/?3^ + (Q^aj+ 

RVC(, q) (Rj_^2 6^Rs 4 ^ Rt + C^Rs 4 /?q 4 ^ ^, 0 4 S^c z) ~ t [ Q, Qj 

4 ^(0 4 qf * s ( qJ + i>^ ') Ri2^ 


R\/(2, l'^ = — .'q j ^3 K) Rq 4 2 (G/bj + Qjb,^ ^10 4 (^QfQj b, b^') (/?i3 + 2 

4 2 f«,’- br ) R,h 4 _ b; ) /?,^ + b^ b^) /? 3 , 4 CgI - b* ) R 32 } 

RV(2,2) =0 

Rv(2, 3)=0 

R V (J, |<) = - A 4 ( i [ R, + 2 ( 0 j i>3 + a3fc,')R„+ CoN b; ) Rj 3 + (a,«,-b, b.){J R 55 + Rj,) 

+ oj 4 ; ) /?„ + fa, O 3 - b, b^) R 3 , + (c| - 4 ) R 33 } 

f^V( 2 ,s)=o 
RV (2,6)= 0 

RV(z,n)=. -N,,-L^(ci\^a2 b,) Rq+2Ci, b, /?,o + 2G^ R„ +Ga,h^R,2 + 2( Q,a^- b, b^) 

4 2(02^3 — b^ Rjq 4 (^,^ 3— b, b^') 4 (Qz^ 3 ~ h. ^s') ^33 J 

RV(2,i)=o 
RV(Zi9)-o 
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+ 4 {(^Z^J 2Ca,Qj+ 3fc,b3^/?,o+(3Q,f^-a,b,)^^^ 

Rv(s,z) = o 

f ^ V ( 3 , 3 ) = 0 

(?Vf3,<,)= -^dUe„ - f^^Z((a,q^^Sk,h,')R^ + 3i,/>3)/?„ 3-30,6, (?„ 63(30, 6,-0,l,V 

Rl/^S.S) =: 0 


f ^ V ( Z , 6 ) = o 

RV(3,7)=: 0Jr„) ~ ^ f (Q,Q, +3b, b,') ('q,"+ 3b^ )R,^ + (q| +3b^) R„ 

+ 3 ffl 3"4 3 i,;;)^,^^ 2 ( 3 Q,b 3 -a 3 b.)R,, 

+ r 2 Q 3 b^-Q^b 3 )R^^J 
RI/(3,8) = o 
RV ( 3 ,‘ l)^o 


Ri^(7,l) = 

RK 7 , 2 ) %^^( a , a ,■^= l ,, b ,) s , i -^( a , Q ^^ ih ,^^) i ,^ 

" *^KK)R,-,-2kK s„-fc,fc3 (S„ - 6 ^,. ^ s^,) _ 3fc,t^ s„- 5^4 

/?l/r76J)= 

<^C 7 •<- <jfr + 5s)| - A/^ X ^3 ('af ^.J>; ^ (S 3 - -t- QiQj-f 2b, b^)(R 3 +S^-Rg) 

«35^ ^<'""'^34-36, b3)S7^r«.«3 + 3b,b3)(R^+Sg)+^?:(S4 + /R4)? 

Rl/(7,i/) =0 ^ 

/? V(l, 5) = - 1 (^ 5,,^ ^ ^ ^ c^ 6 ^, R,^) - % i(- bf S„ - b, b 3 R,, )S„ 

+ C«.«3 ^ 3 b, b3)S ,, 03 ^ 2 b^ b^) /?,^ - b, b^ f + 5,,) - bj /? ? 

1 / 7,^) _ + 4 Rz + ^Yr^ + Sg -|o) + (S 4 + /^) + 3<^V/?v - i^) ^ 61^ i-miRz+S^)] 

4-S,) + ('2a,Q^^.3b,b,)(S4 + R,)+3(Q^b^)(^(- R^) + Q^Rs 

RV(7,7)=zo 

Rv(7, %) = -±{p^ - (f 

+ +3br)S,„ f r«,C7, + 3b,b,)(f?,,+ S;„) + fq,Q3 ^;j/,,b3) + 

^ ^ ^ ^ + ‘PA^S f S 7 + g <^,(^7 +Sg)} - AiX 4 {(2 a, Q3 + 3 b, b 3 ) S 5 

■*• ( 20 ^ 03 + 4.q2s^ + ^ 
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RV(S,1) = - + 

+ 2 («| 63 +Q 3 b,)S^H- ^('Q,b 3 + a3fc2X/?T + S^.) + -^(^i^ii-b,t>j)(S,4+ 2 S,£') + (^a3-t^)S2i • 

+ ^(a,a 3 -b,b 3 Xs ,7 + 2 SjoV<|'^ 2 ^i“^=»*' 3 ')(^'‘' 

R\j(<i,z) = 0 
RV(%i3) =. 0 

RV( 2 ,^) = - N^{ ^Q,b,(^ 3-|2 ^■s:^) +la,h^(R^-^) -^<^ 3^5 

+ ^ K t>3 + Q 3 b,)(R^ 4 - Sg) + 7} ('Q<^- fc.") S,4 X(Q,a3 ^- b,b3 ) (R,t + V S2,) 

Ri/re.s) = o 

R\/(^,(,) = o 

/? V(9,7) = - W^fa «'b,S7 4- f fQ,^3 + c23b,)Ss4^a2b3+ Q3b^)/?s + ll-(atb^+ +S&) 

+ V rQ/Qi- b,b,)(»?,7 4-S,7) %^Qx-b2)R,q^-it<'«.«3-b,b3)S,2-^ ^ 

■^^K-b!)Sa} 

RV(9,^) - 0 

RV('s,<i) =0 


Rl/^q, ( ) = -{So 4 f(S^~ |a) -<-^<;b (/?3- & +Sj.) + ^(52, 4-<^6{.r^74-Sg)| - fVA ^ 

4- a(a.a, +3b, bJ(/?3 - ^ + s^) + («| + 3b^ )(s<, + /?6) -^(al-i-ibl) >2(a,a3+ 3i>, b^)Sj 
4 ( 0^03 + ^b^k')(R7^Sg)^(SR,h^~0^b,)S,^ + (3Q,b^-a^b,)S,^ 4 2f3ct^b -d, S,g 
4 (3Q;jb3-a3b,)S,q + 2(3Q3b,-a,b3')S:^^4(3t?3b,-aAX^ao^^30-^^^3b3S2i} 


/?V(^2') = 0 ^ 

4 3(a^4 3b^ )(/?4 - la) +r<?3 3 bJ)% 4-(a,Q3 4 3b,b3)(R^ 4 Sg) 4 2Q, b, s,^ 4 3 (a, bj- Q3I),) R,^ 

4 (ZQ^b^-a3b^)R,cft-(3Q3b,-^tb^)(R20-^^zl)-^^^i^3^i^} 


RV(‘1,5)~0 


Rl/(' 9 i 6 ) =0 

RV(R,1)=- ('^ ©, 4<5b©^ Rg) _ Nfl . ^(2fd,Q^ 7 - 3 b^bj^Sg 4. 2 f + 3b.b3) R^ 4(a,V3b')S^ 

4 4 3b,b^)(/?7 4 Sg) ^ z a, b,S,j 4 f ^Q,b^-Rzb,')R,^ +( 3 a^b, -Q.,b^) S,c( 4 2Q2b2R,t, 

4 (3Q.-^bf — R/ba) 5^2 4 (3Q.jb^-- ^ 2 ^a) 


R\/(R,%) = o 
RV(^>R)=o 
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RVO(0=-o 

RV0(2) = =-^/iv[2Q,b,/?3 ■»-{Q.i’z + a.M^fe + (Q,^4-a3b,)R^ 

+ («,t3 + a3b,)/?g + (a, a^_!,,t. ^ r,^) + (aAo,~kh^)iRn + R^o') +(al-t |)«22 

RVO(l) --(Ro-'r ^5 +/i<^ f?t+ |!S ^/?7 

+ (G^ +rQl+ 3b3')R54-(a,a, + 2 />,b^')(?^+ ('a,a 3 + sb,bj ))?7 + (3a,b^-a^l^)/?,^ 

+ (3Q,b3-Q3b.)f?,^ ^('3a2b,-a,b,)R,g +(3a^b3-a3b,)R,c, +(3a3b,^a»b3Xo+«3^j^a^ 

RVO(s)~ zQ,b, ^3 + + ctQ^b^S^ +(a,h^+azb,')s^+(a,b^+a.3l>i)Srj 

+ (®a^3 + *^ 3 ^ 3 ) % ^1 K") ■*• ^Ig) (^1^3 ~ ^3)C^(7 ■*■ 

h C^S-bl)S22 -f- ~ kz^sY^'^f 

R I/O a ) = + ^"$3 + + 4" S 5 ^ ^ I 6L5^ +<f, (Pc Ss A/;, -i { (a,=^ + 3b? )S 3 ■*■ (al+ s )S., 

^ +•3^3* ^Ss + (d,Q^ + Sb.b^Yi ■+ r^.G3 +• 3b,b3)S7+faj(l3^3febj)S3^-f3Q,b^-a^^^S'^+(30,^-a3b,)^ 

+ f3 G J b, - rt, b,) S,g + (3 b 3 - ^ 3 b,) S,«, + fi « 3 bl -a, b 3 ) S;J 0 4 (3 tt 3 b^- Q, b 3 ) S^, +. ? Q 3 b 3 S«} 
RVO(l) - 0 

RVO(S)= i;{=iR,b,S, ■>- (<i,b^-^Q3k)S^■^{^^b^-l■0sK)Rz■^ (Q.,^‘-bf)S^^ -v C^^.^z-b^ b^')S,j, 
4(Q,G3-b,^)S,5 4 

RV 0 (‘?)= 4 <ib(9^/?^)-./v^i[(aN5b;^)s, -f r«.Q 3 + ^b,b3)s^4 C«^aj+3b,b3)<?z 

+ 2 C,b, 5)3 -Q,b,)s,^ + P«3b,-C<b3)S,s t^Qib^ /?,<,} 


Note : 

RV(i,j) for i and j varying from 4 to 6 can be obtained from 
RV(i-3/ j-3) by replacing Rj^ with -Sj^^ in the expressions 


\ 
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Periodic aerodynamic forces terms for Parametric Resonance 


i) 3-^. i + 4. ta.b, R;^^ 

+ 7jfQ,b, + Qab,)f?je + 

Rv (l, 2 )= 0 

RV(i,i) = 0 

l^\/ (l,u) ~ ^ 3 ) ~ i ^-I ^'1 ■'■ ^ '^-^3 -‘- ^ r®' ') S 

4- + a,b,')R ,7 4- 4- 4'(^I^3 +^3^, )^3I •i^‘^3^33 ^ 33 } 


C'/S) = 0 

/ei/ = 0 


R\/(l, 7 ) = 




RV(I. ?:) = o 
RV { 1 ,^ 1 )= 0 


RV( 2 ,i) - o 

/? V C5 , 3 ) =: - ^ f R>, 3 + '?«6 ^(7 ^/! <^C i ^4*4 €^^25 + <S?/ “ '^A { «< ^5 3 

4-(Z| /?,s +^i 

RV(2,^) = o 

Pi/ (2,5)= 4-i(9,Pg 4- if (9c 4-<^(9,P„ / 

R V’ (-?; 6 ) = ~ h{^>H ^ ^'S ^(«t 4- 4- Pjp) — ^<3 (Qj Rj 7 

(^:?g 4- 4 C 3 P 30 I 

RV (3| 7 ) = o 

RV (' 2 . 9 ') = + ?(9c ^ ^7 ^ ^'R.,4-3®; ^, 2 ) 4 - A/^ -L (a,^ P.o 

+ Qz P„ 4- Qj P,2 } 

» 

f?V(2,q)=-ifP,s-^/5P,^+C^ (9cRji + P.<f /?3.4^c«c'?32+‘fe®c^33)-/VA-Ji{^.Ga'^3, 

“*■ ^1^3 ^32 + (^ 2*^3 ^33^ 


I 
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RV(^,^) = o 


/?1/ f ^ ^ f R,J 4 PA, /?,3 4 1 Rx^ "■ - Na U . 

4 bi Rai^ ■+ bj Rj 5 + bj ^at,| 

RV(3,3) ^ 

RV(S,M)= o 

R V (S . 5) = -1 (R.^ + a ^ ^ + i Rn + <Sc ~ '^A M 

^,9 - Rac, 4 b; r4 

RV (3, 6) = - ^(acja R4 + H R^ 4 <9, 4 Rci + A<fe(9c R.^ 4 ^0l R,a') “ ^a{ L R,a] 

RV(i,7) = o 

(3,%) = -r; ( '^‘S 1^ Rxo-^ t Ri2L+ §^^R3l-4 %S<,R35.4<1 ^^®c Rsa) - '^A-q{b,bj(?3, 

4- b, bj ^32 4 Rjgl 

RV(^,R) = - ^fHj.4 3S,(?s 4 -^^Rt 


Rvr7,0 = f^Aif{^(^<-b")(h-^)-^z(<i,Q^-b,bA(Rs-^-f-S^)-i-(al-hl)(s^^R(,) 

+ bl )S^^- ?f«/Q3- b, b^)S^ -f- CQaas- ba b^^Ry'l-S^') - ta.4Qab,)(S, +^ 

- f C«/b 3 f 03 b, X^,7 + 5 Sao)-|(Qab34- a^ bj)('S,q + R,o 4S;,,) - Q^b^Saa } 

RV/C7;2) = 0 
RV(7^ 3) = 0 

^ 1/ (7, tv) = /Va ^{(qA b; )(R3 - ^ 4- S^ ) + :i('^, Qa- h, h^)(S^ + R^) 4 s(q^_ b5(4 - ^) 

4 (Ri-bl)R^ + r«,Q5-b,b3)(R,4Sg')- a,b,S,^-i(a.b34a3b,)(R„+R3o4Sa,) 
"iC‘^ab3+ Q.3bj^)Rjq —^abjRaal 
RV(7,5) = o 
RVO>(>)~0 

Rv(i,i)=. (q ,«3 _ b, b 3 ) S 3 ^ ^(a,a,,~bAs) ^S + CR"- bf ) S 7 4- fi.^x-b,bAR;^^) 

-«ib,S ',7 - 5(^1 b^4 aab,)(R,^ 4 ^iq') —^aba R,c( - i(^ib3-v<3.3b,')S3j 

— Ka^bj-v 0.3 b^") Raa^ 


RV(7,s) = o 
RV(7,R) = o 
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/ A-0 <- N^.3Bct3(Rs-% % Q Q3 (Ri^^s')J 

»>“'-;i;;cA*.-‘.«*t'':T-t!-:-.-*‘'-‘'’"‘‘-"’ 

I J S S,3 + e t '& 

“"’ it,.. ..«..«<.'i*‘' i.;-«>v--.'-' . , 

\ o o 4.6^5- 4- 

Rv(«n^=^^ . s ' 1 , I la!s„ 

« «-!- ‘ ' 5 Ntpfl, Saints. 


RV(^> 0 " 

RVO< ' 
RVC^/^) 


+ 61^5 , 

^|(S),(R7+-^^1' ^ .2^ +b,h,CV^^ 

=^0 /r? xS ^ + ’ 


rvC^.s) 

KVO.^'i 

RV(^.'^V 
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RVO(i)= + )^5 +(«I«2.-^(4)/?6 

-■^('«/b3 4a3b,)<?;jo ~ ■i(‘^3.i>3+^3K)R:il — Q3b3^Zl| 

RV0(2)= 0 

KV0(3) = o 

RVO(h )='NAi ^{ ^ 4 (Q^-h^^Ss -y,bjs^ +(«.Q3-i>.0^7 

'*' ^ 3 ) ^s- ■'■3 ^3 ■*■ ^Z^i)^i(, " "a (‘^ 1 ^+ ^n" "a 

“ -I- - i (^ib^-*-^3^i) ^zo~i(^i^3 ■f-^3k)^zi --^3^3^zz } 

RVO(S)=0 
RVO (6) ~ 0 

Rvo(n) = NA-^{(<i^-b^)s, +(R,Q3-b,h^')S^+(a^Q^~b^h^)R^ -i-C«at>,4 

— b, S^3 — ^ (cii s,^ — ^z ^z ^'‘1 1 . 

RVO(S) == O 
/?l/0('‘)) = 0 

Note ; 

RV(i,j) for i and j varying from 4 to 6 can be obtained from 
RV(i-3, j-3) by replacing R^^ with -S^ in the expressions 


\ 
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APPENDIX VIII 


The elements of matrix AA for flapping-lagging-feathering 
rotor consist of tv;o separate parts 

(a) contribution from inertia and stiffness, AI 

(b) contribution from aerodynamic forces, AA 

Thus 

A^j = AI(i,j) + ^ AA(i,j) 

A^q = AlO(i) + ^ AAO(i) 

The various elements are given below; 


— l-t-e — 6^- j + y! ) ■>- + ;^fi + oc’)(Q 3 4-b*'^ 4- 

= 0 

flr(l, 7 )= 
fli 0 


91 ( 3 , z)- e-N^{i(a,^+b^) 

Pil (3,i)= f^j{oca,Q^j 

91(3, h) = o 




9i(3,s)~ od^a, 

9 i( 3 ,q~o 

9I(3,'7)~ ol^ 

91(3,^)- + Nj_{-Lcc(3a(+ t*) ^ 
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+ h(^3 ^*=" 3 ')! 

AI(3,4)=-0 
fil(3,5) := 

flr( 3 , 6 ) = 2^2 

flr(3, 7 ) = O , ,2 I 

f]I (3,‘i) = ^JI-CC^)-^^ ^c " '^l{i('~^ ^ 


Aia2) = 0 

flr('f,3) = o 2 ,2, 

Ar(<i.s)= 
fli(^,6)= o 

AI0,‘?)=^ 


+ ^ 3 )} 


m^>o=° , j -■) 

Ai(s-,2)= '^li 


fir (s,‘<)= 0 

fll{'5,5)= 
fiJ(S,6) =0 


fi r (Si'j') — ® 
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PI(t,0 = 0 

qi(i,ii) = o 

flr(6,s)=o 

R I ((>>(>)= zo-z^ 

fir(^,^) = 0 

fl I ((> 1 %) = - ^oc(R,^i9^ -I- ~ ^ cc a,a 2 ^ 

RI((>iR)= 

fll(7,Z)= 0 
Rl(7,2) = 0 

/) r C?, 4) = A/j [oc + 4?2 Q^l 

f\l(7,S)^ HV(9-<5s')^m 
flr(7,6) = iNj{aii)jR»w} 

^ ^ ^7, 7 ) = /?^ - /S % Nj. { 1 ± ( af + ) j 

ftl(l,<i) = 0 
flI(l,R)^o 

fll(i,\) = (/+oc*)«i 

RJ(8,2)=^ 0-oc^)<p^-2^cp^+N^^jb,b^(K-0] 
f\I(8,2)=^-Zcc^^^N,{{pc(il- -^oC(a,^4. b^)j 
ni(i,R)= 0 
PI(8,5)= w//?^Q3 
/)!('?, 6)= 
m (2,y)= a<jb^?2 

Air r?/R) = ^n, fz-oc^) (l-c^)4-N,(^0 + c^)(3a,^^l>^) + 
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= 0 

=■ 2CC ^ -I- ^ (0^+ b, ) + j 

/}:n,3)= r/-c^*)<^ + C^’_,)q,q3| 

Pj(q,q)=0 

fli s') — COp 

flr(^0= =^Pct^z 

f\l (^,n)= HocRr^a^ 

fii C9,9)= R^(,q 

ft I (^/‘^) = Rrn(\- OC^) — (3 (l -h ec^) + <^ ^(_ OC ) + |"- -^(*/ + 3 b, )('l + eC^) -4- -|^ (l- 3ec ) 


fliOfO = N,{(i^ s^'")^za,] 
RI0(2) = C 
flrofs) = 0 

f)IO(H)= Nj{^Ci-cc*)(a,Q3 + l^b^)} 
fliO(s)= (e-<K^-^\ 6 ^^)a 2 


RIO(&) = 0 

RIO(l) = +/^r{(2~ RryiCC^)a,Q2} 

R[0(8)=0 


RW(R) = o 
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Aerodynamic contribution terms 


A A f I , I ) - - (9, 1 a, 63 fa, a 3 4 1 , b,) n, b, + oc o, b, - 1 F3, 

+ Foe fcJifij + b, Fjg •*■ T F j-Aj 

fl/^(l, 2 )- 0 

f)R{h^) = ° 

- - Wa f ;ffa,P3 4 b,b^) F, 4 0,^3 F„ ~oca,b, ^ pcq, b, F„ 4 ^ocfajb.-Q.^^p^^ 

4 F ^ 2^3 ^7 ~ bj 0 - 1 Q 2 bi| 5 I 

ftf\(t,S)~ 0 

/?/?(/, 6 )=. - Nf^f--La^b^ (F^ +±fq,b 3 _ + 4 ^ab, F,,g 

- F a^ b^ I 

^/9f,,v)= -(F^4^’/^, ^-{KFb;^)fio+i:«2 f^, + |K + b])ha 

., orfO, b,- Q.fc, )F,^ tj + 2 (0 jl, _ S,b3^ Fj, - f a.bj F3,} 

/)fl(l,s) = 0 

ftf)(i,‘i}=o 


f) fq , I ) =. — ('cZj ^ 4 (^a 2 ^ + ^'i ^*^2 ^s') 

f2i 2) = <■ ’? 5 * ^ + <a fi * ^ ^€<? fio ■► 5 SV,s)- « a{ f '"1 * ib,'’. <=,o 

+ ^ ec b, bj ^ ec^QiQj ■!■ X 

/)/)(-J,3) = K(fV3 *<$ »«< (7, + ,^3 * l’<iV b «■’ F«) <- ^aH “■‘’j'? * f “'"a '53 

+ f ro,Vfc,*)i 5 , + A«ca.’ 

nfi(i,'i)=-ta,(F^*s,F„) 
fift ( 2 , 5 ) = 0 
flR fq,6) = 0 
Rf}(^,7)=F a,F^ 

4 ,f 4 .^F, 4 ^<j 6 ^ 

+ 3a\F,,+i(al + bl ) + c^Ynf- b^ ) F3^ 4 F3^ + a^ b, F3, J 

RR (^1^) - cefFs ^^ 2 o-*-%’iir''^cFz+ r’^<f^F 2 ,-*- ^ < 9 cFz 2 + ^ F‘i)~ '^A{f ^1^3 <3203 'so 

- f :ca.Q2F3, - f a,c?3 ^ a2«3 ^3 ~ f F,} 
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f}fl(3,l ) = - 


nn(3,z) = ^ 4f^s^^cF,^) ~ ^//i. { ) 'iv 

+ 2 (a,\- 36l)F,,^f a,’ F„ . f (al- fc/ ) F,, ^ ^ «,Q, h,3 + ^ q.q, F, J 

fl fl fs, 3) = - (F; + ? ^ F 3 + (^ Fg + (9, 4- ckfi^ F, + :?P^ g. 4- C5(>V' ^) _ /\/4 i Q,Q3 F, ^ 


flflf3,q)=o 

Rn(3,S) = o 

Bfifs.O = - a^( /;, + 1 sj (J, ♦ /r„ + si/ + <j’ ij, ) - ^ «/ f, | 

f\f]{ 3 ,n)=. cca^F,^ 


Rn(3,S)= -oc(F,^ + qr^ + ^,+^<9c F3? '■33 4-c^"F 3<,)- WA{-f b,t>3F^6 

' •*• f f ^.^2 'Si + f tj 'Sj. "• f ^ 2^3 'S.i -f ^2 ^S^i } 

FF (3,H) =: ~ (F^ ^ Z Q^F^ -i- <^ Fg + + l?<^ f * fS> + -F 3(9c f);^) - F/a { -q- Q| S 

+ ^ 3i>;) fr + X q/ , 5 ; 4- a « a,i >3 F 3 , _ ^ a,b 3 F 33 

- fa, ^,F,^-f a^h, F^, - ^^oc^o-l I 


flR(‘l,>) = <^, f-<ib(?^Gc, h-'Va{i‘^ 2«3^^ +(«.a3 +b,b3)G,o-f 0;jb,Q^32-oca,b,(5;,g-«a,b^G3, 
-f C?3 -f b, b3) G3^ * i ^3 Gs,s, I 

^/l('i/,2)=0 

flflfV;3)=0 

flfl(q,.‘/)= F (C{,<3.^ + b,b^)Gc^ -f- a^Qs^o -«^<^2^, ^25 +f 

+ +fV«.F3 + b,b3-)Gq3} 

flfl(i<,5)= O 


fiR(ti, 6 )= )V^ [-{^z^.Fik, -f (<3,63-03 b,')G 17 + i<?2 6, Gig + i^^bj G,q + i(^^h,-aib^')G,:,o 

~ i j 


f\P(H,^)= 4 ^*G,o 4 . 


<^*G„4 6^*G,z + ^A{i(^|Fx')'^q +-LK+b" 'JG.O +^«xG„ 

4- CCC0,b3-Q3b,)G^^ +^«,b3G3^+ |(fit3b,-a,b3')G32- f-«a^G3iJ 


f1fl(^,£)= O 

ftfl (^y,9)= o 
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+ j- + -J- rc d2.^j(3-43 j ( ^ „ r, r 

/lfl(5,j-) = -‘<(S/3*1i‘^«*'‘‘”'i ^‘- ” " , ^ . _Sa,a,G„ 

3 S (o; - 5 G3j--il “»‘>3&-,G- 3 ”="“3^^ 

ftAC^iM)— X<3a(Gq+ <9c G|() 
flflCS.S) = o 
flfl (S-,6)= 0 

/)/) ('5,7)= 0,^3 ri, 2 i^^r 

/}flrs.«)= ' 

4 . i- «r 5-, G,. - 4 oc ci,<^2 G 3 , — ^ aiQ, Gji - -£ f 


Z','i = <C i «» * ^ f <=- : ® 

4 K{a,*4.fc,")Gj^ 3^ff,'Gjs3^(a| +3fc3)G;:4+fG4 S„* ^ . i 3S 3 

3ic<G,t,G,, 3i«>«,«3<^36 4lA»3<;„-5’‘’3“3S-3^ 

H 

ftflC6t^) = ° 
flf\(^,s)=o 

rtfl (^, 6 ) = Qz (G,^ -^n^zi -^ ■" V ';'a { 4 ^7 

ftfl('6,g)=OcCG,s + P.G2o+tcG„ + ®cG« + l‘^^3( ^-^c®c^3Z + ‘fe®c<^33VWA{ ^ 

f (G3o + G35>ja3-^^^'^2^31 b,bj&Sl} 

fiA (6,q) = <^J+ 5 (S(:G5+ ^ G 7 +«^Gg +g<^Gq +f |J|0 Gi,| + 3IS5: G,2 

^ S -.f (a| + b^)G,,~ f <^ab, Gs, -4- ^30- "^ 53 ) 

4jY‘^36 ^57' 
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W7,0 = - Go- «. ( E. <- G,) - <-S) - I Y-Si - Ss ) - ff/Gs - »!, I i 4 G,') 

+ (ci,a, -f b, b^)Gy +a,a^{F.^G,-.^)^^al (f, 4 | i>’ ) (g^ - 

+ ^Gs-fa^b^G,^ + f r«.^3-‘^3^)G,7 ^^Qzb,G,g -+ f Q^i>3G,c, - |pa,i,3 (Fio + Gy,) 

4 ^Yq, aj 4 fc, g,/2 4- } 

(]R(1,2) — O 

Rf)(l,z) — o 

f\R(l,L{)= ~F^ - /^A(ifa,Q3-+ h,h^)(F^^Gg) + a^Q^ fT 

. ^ ^ ^ I «.V - f ) 4 i ra3% ^3^ ^ 

~ fi, f 4,q, F3^ 4 F3^ /^gj 

f\f\(-i,s) —o 

M(i,i.) = f _ S,^-) , -±B,I>^ f;^] 

flfl(7,7) = -g G^ ia,a,(r^t<Sg) j-J('a,7+t,’-)G, + 

4 (^< <^3 -I- b,i^)G^ -i- f(~ Q^b,) Fly ^ ^Q^b, 1 

flflC7,S)=.o ^ ^ 2 ^ 3 

Ffl(-r,s)=0 


RR (Sjl) a^f- 4. (f^ 4 Gs) -*- ^2 f ( *^5 ■+ ^ 4 ^ <^2 'fe ( *^6 + G 

flfl (S,2 ) =-Go- 3|G,-<^F. 4(9cGi-«^(gc(S4G^')-2^cScG7-3^<^CF34G6-S)-c^*(Fg+GH^ 

- -'Va{- 2 jQzQ 5 (F, + Gs) 4 l„Q,a^(F34G^- ^')4 Q b 3 &^ 

+ ^°'(F,4G4^ 4J«4bY(G3-|o)-4-li(3Ql+t;')Gs4^YfF3^4i*(ar4bYG,6 

+sYzG 3^+ ^V^5 + 3l?3)G4o 4 & b,b3G4^| 
flfl(S,3) _ ‘^('‘feFis 4^ G,3 4 ^c^G,^ 4^ ■^'fe<^cP.7 4|(S,&n4g*Gm') - 

- ^ («^ bn G„ + it. ( p, 4 GO * z «.(,, G, - s ) _ i ^ „ a, a, f;,- f ».t,S 

f a, «3 G„ - a (T^ 4 = a, a. G„ 4 ^ a, G„ 4 a („j . bj 5 g ,,4 ? a.c,, G„ - «V I, a t 

4 3C,<;bC^'_ 

flflCSiS) =0 ^ 

/?/? ("g, ^ ) = o 

(\R(2,i) = 4Cjfii. fg 

flfl (S,S)=-{i <5a + 4 /3<^ ( F, 4 Gs ) 4g*(S^ 4<^^ 4 3c^(? 4 3;8(9^ Gj }- ^ (sar^ bf •) G^ 4 3 q^ 

4lb,b3Gs 4^a,b.(F3,4G^,) F3^_4^Yb3C;4o4^\a,^_b^)G4^ '' 






r Awl « 

«„r.,e)=-"(^‘^-y‘ " ’ %, .Vr.G,^ 

4 Ko’ R6'^Tr“'‘‘* . ,fa&i+s|'(j('i-*‘^‘''?^''^‘ ‘ V 

"m 

flfl(^,H') = -ocQiftM 

fifi('i.'')= ° , a + «(a,V3b;-)a„ 

+ iS(a, -b, ) n M ■) J j-<i,aj(bT + ^s) 

Ana'i') = 't'^ ,„2f 4 - 4 <ja 5 +ia'Fi35s H»' 

j.J-fq?t3b’')A + 5‘‘= F* , _S.o,KF,i\ 

h' ‘ , . oClVo,-t,)Gi,2 H. 2 3 1 i 


£C 

H 


"1 

S „ /Xra,'*bf)F3*2“*'^ 

,3 fF + g=fi*'S ® 4 » F 4^0,0.%.*- 

4,.G 

RP0(^)-° ^ 7 ^ ^ o Gf,, vQz°'3^iS 

,„o (. ) -. ^ € A ^ S= 4;, 4 . t ^ f 

4xc«;4i4',«s*^;4 ‘ '<^5 
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FORCED SOLUTIONS FOR CASE 
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FORCED SOLUTIONS FOR CASE II 
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FORCED SOLUTIONS FOR CASE III 
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FIG. 3a FORCED OSCILLATIONS, CASE HI 
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FIG. 4 PARAMETRIC RESONANCE, CASE I 
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FIG. 5a PARAMETRIC RESONANCE, CASE H 









FIG. 6a PARAMETRIC RESONANCE, CASE HL 
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FIG. 12 PARAMETRIC RESONANCE WITH AERODYNAMIC FORCES. CASE I 
























FIG, 18 FLUTTER SOLUTION (EFFECT OF LOCK NUMBER Y) 
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FIG. 20 FLUTTER SOLUTION (PENETRATION OF INSTABILITY BOUNDARIES 
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FIG. 23 FORCED OSCILLATIONS OF FLAPPING-LAGGING-FEATHERIMG ROTOR 
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